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Abstract 

The mclusive e'^p single and double differential cross sections for neutral and charged 
current processes are measured with the HI detector at HERA. The data were taken in 1999 
and 2000 at a centre-of-mass energy of ^/s = 319 GeV and correspond to an integrated 
luminosity of 65.2 pb~^. The cross sections are measured in the range of four-momentum 
transfer squared between 100 and SOOOOGeV^ and Bjorken x between 0.0013 and 
0.65. The neutral current analysis for the new e^p data and the earlier e~p data taken in 
1998 and 1999 is extended to small energies of the scattered electron and therefore to higher 
values of inelasticity y, allowing a determination of the longitudinal structure function Fl 
at high (110 — 700 GeV^). A new measurement of the structure function xF^ is obtained 
using the new e+p and previously published e^p neutral current cross section data at high 
Q^. These data together with HI low precision data are further used to perform new 
next-to-leading order QCD analyses in the framework of the Standard Model to extract 
flavour separated parton distributions in the proton. 
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1 Introduction 



In 1992 the HERA accelerator began operation, colliding lepton and proton beams within the HI 
and ZEUS experiments. The phase space covered by HERA inclusive deep inelastic scattering 
(DIS) cross section measurements ranges from small Bjorken x at low Q"^, the four-momentum 
transfer squared, to large x at Q'^ values larger than the squared masses of the W and Z gauge 
bosons. These measurements provide an insight into the partonic structure of matter and the dy- 
namics of strong interactions and test quantum chromodynamics (QCD) over a huge kinematic 
range. 

Both neutral current (NC) interactions, ep — > eX via 7 or Z° exchange, and charged current 
(CC) interactions, ep uX via W exchange, can be observed at HERA, yielding complemen- 
tary information on the QCD and electroweak (EW) parts of the Standard Model. The cross 
sections are defined in terms of three kinematic variables Q'^, x and y, where y quantifies the 
inelasticity of the interaction. The kinematic variables are related via = sxy, where s is the 
ep centre-of-mass energy squared. 

Measurements of the NC and CC cross sections in e+p scattering have been made by HI and 
ZEUS based on ~ 40pb~^ data sets taken between 1994 and 1997 [1,2] with protons of energy 
820 GcV and positrons of energy 27.6 GeV, leading to a centre-of-mass energy ^/s = 301 GeV. 
Here, new e^p NC and CC cross section measurements, based on data taken at y/s = 319 GeV 
in 1999 and 2000, are presented with improved precision using a luminosity of 65.2 pb~^. 
The increased centre-of-mass energy stems from the change in the proton beam energy from 
820 GeV to 920 GeV since 1998. These data sets together provide the most accurate neutral 
and charged current cross sections measured by HI at high (> 100 GeV^) in the first phase 
of HERA operation (HERA-I). 

The NC analysis is extended to higher y up to 0.9 for 100 GeV^ < < 800GeVl This 
extension of the kinematic range allows a determination of the longitudinal structure function, 
Fl{x, Q"^), to be made at high for the first time. This analysis is performed on both the 
99 — 00 e+p data and the e^p data, taken in 1998 and 1999 with a luminosity of 16.4 pb~^ at 
= 319 GeV. The extended high-y e~p analysis and Fl extraction complement the inclusive 
cross section measurements published in [3]. The difference in NC cross sections between e+p 
and e~p scattering at high Q"^ is employed to update the measurement of xF^, superseding the 
earlier measurement [3]. 

The accuracy and kinematic coverage of the HI neutral and charged current cross section data 
enable dedicated QCD analyses, which test the prediction of logarithmic scaling violations over 
four orders of magnitude in Q'^ and allow parton distribution functions (PDFs) of the proton to 
be deduced. This in turn allows predictions to be made for future facilities such as the LHC, 
deviations from which may be due to exotic phenomena beyond the Standard Model. 

A next- to-leading order (NLO) QCD analysis of the HI data alone is performed, using a novel 
decomposition of the quark species into the up- and down-type quark distributions to which 
the NC and CC cross section data are sensitive. The fit parameter space is narrowed using 
theoretical constraints adapted to the new ansatz and the experimental and phenomenological 
uncertainties are systematically approached. This leads to a description of the complete set of 
NC and CC data as well as to new determinations of the PDFs and their uncertainties. For 
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comparison, the QCD analysis is furi;her extended to include the accurate proton and deuteron 
data from the BCDMS muon scattering experiment [4]. 

This paper is organised as follows. In section 2 the definitions of the inclusive NC and CC 
cross sections are given. In section 3 the detector, simulation and measurement procedures are 
described. The QCD analysis method is explained in section 4, followed by the measurements 
and the QCD analysis results in section 5. The paper is summarised in section 6. 



2 Neutral and Charged Current Cross Sections 

After correction for QED radiative effects, the measured NC cross section for the process 
e^p — > e^X with unpolarised beams is given by 



with (t)%c = Y+F^TY-xFs-y'F,, (2) 

where a = a{Q'^ = 0) is the fine structure constant. The A^'^^^;"^"'^ corrections are defined 
in [5], with a and the Z and W boson masses (taken here as in [3] to be Mz = 91.187 GeV 
and Mw — 80.41 GeV) as the main electroweak inputs. The weak corrections are typically 
less than 1% and never more than 3%. The NC structure function term 0^^^ was introduced 
in [1] and is expressed in terms of the generalised structure functions F2, xF^ and F^. The 
helicity dependences of the electroweak interaction are contained in Kt = 1 ± (1 — y)^. The 
generalised structure functions F2 and xF^ can be further decomposed as [6] 



^^3 = - ae ,J^\,2^ ^Ff + (2^eae) ^F,' , (4) 



(Q^ + Mir ' ^^^^^^ \Q^ + M: 

with K~ = 4-^(1 — -j^) in the on-mass-shell scheme [7]. The quantities Ve and are the 

vector and axial- vector weak couplings of the electron^ to the Z° [7]. The electromagnetic 
structure function F2 originates from photon exchange only. The functions F2 and xF^ are 
the contributions to F2 and xF^ from Z° exchange and the functions and xF^^ are the 
contributions from interference. The longitudinal structure function F^ may be decomposed 
in a manner similar to F2. Its contribution is significant only at high y. 

Over most of the kinematic domain at HERA the dominant contribution to the cross section 
comes from pure photon exchange via F2. The contributions due to boson exchange only 
become important at large values of Q^. For longitudinally unpolarised lepton beams the F2 
contribution is the same for and for e+ scattering, while the xF^ contribution changes sign 
as can be seen in eq. 2. 



^In this paper "electron" refers generically to both electrons and positrons. Where distinction is required the 
symbols e+ and e~ are used. 
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In the quark parton model (QPM) the structure functions F2, F^^ and are related to the sum 
of the quark and anti-quark momentum distributions, xq{x, Q^) and xq{x, Q"^), 

[F2, Ff, F#] = X 2e,v,, v\ + ajjg + q} (5) 

1 

and the structure functions xF'l'^ and xF^ to their difference, which determines the valence 
quark distributions, xq^ix^ Q'^), 

[xF^^, xF^] = 2x ^[e^Og, Vqaq\{q -q}^2x^ [e^Og, Vqag\q^ . (6) 

q q=u,d 

In eqs. 5 and 6, is the electric charge of quark q and Vq and are respectively the vector 
and axial- vector weak coupling constants of the quarks to the In the QPM the longitudinal 
structure function = 0. 

For CC interactions the measured unpolarised ep scattering cross section corrected for QED 
radiative effects may be expressed as 



dx dg2 27rx 



.x-lQ^ + M^J -^c(1 + A--), (7) 
with (P^c = ^{Y+WtTY.xWi^-y^W^), (8) 

where Ap^'^"''' represents the CC weak radiative corrections. In this analysis is defined [8] 
using the weak boson masses and is in very good agreement with Gp determined from the 
measurement of the muon lifetime [7]. The CC structure function term [1] is expressed in 
terms of the CC structure functions W^, and xW^, defined in a similar manner to the NC 
structure functions [8]. In the QPM (where = 0), they may be interpreted as lepton beam- 
charge dependent sums and differences of quark and anti-quark distributions and are given for 
an unpolarised lepton beam by 

W+^x{U + D), xW^ = x{D-U), W^=x{U + D), xW^^ = x{U -D) . (9) 

Below the b quark mass threshold, xU, xD, xU and xD are defined respectively as the sum of 
up-type, of down-type and of their anti-quark-type distributions 

xU — x{u + c) 

xU — x{u + c) 

xD — x{d + s) 

xD = x(d + s). (10) 



For the presentation of the subsequent measurements it is convenient to define the NC and CC 
"reduced cross sections" as 



&NC{X,Q ) = — 



1 Q'^ X dVivc 



Y+ 2na^ dxdQ^ 



27rx 



dxdQ^ 



(11) 
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3 Experimental Technique 



3.1 HI Apparatus and Trigger 

The HI co-ordinate system is defined such that the positive z axis is in the direction of the 
outgoing proton beam (forward direction). The polar angle Q is then defined with respect to this 
axis. A full description of the HI detector can be found in [9-11]. The detector components 
most relevant to this analysis are the Liquid Argon (LAr) calorimeter, which measures the posi- 
tions and energies of particles over the range 4° < ^ < 154°, a lead-fibre calorimeter (SPACAL) 
covering the range 153° < 9 < 177°, the Plug calorimeter covering the range 0.7° < 9 < 3.3° 
and the inner tracking detectors, which measure the angles and momenta of charged particles 
over the range 7° < 9 < 165°. In the central region, 25°<6'<165°, the central jet chamber 
(CJC) measures charged track trajectories in the (r, 0) plane and is supplemented by two z drift 
chambers to improve the 9 measurement of reconstructed tracks. The forward tracking detector, 
9^30°, is used to determine the vertex position of events when no reconstructed CJC track is 
found. 

The ep luminosity is determined by measuring the QED bremsstrahlung (ep — > ep^y) event rate 
by tagging the photon in a photon detector located at z = — 103 m. An electron tagger is placed 
at z = —33 m adjacent to the beam-pipe. It is used to check the luminosity measurement and to 
provide information on ep — > eX events at very low Q"^ (photoproduction) where the electron 
scatters through a small angle (n — 9 < 5 mrad). 

NC events are triggered mainly using information from the LAr calorimeter. The calorime- 
ter has a finely segmented geometry allowing the trigger to select localised energy deposits in 
the electromagnetic section of the calorimeter. For electrons with energy above 11 GeV this is 
100% efficient as determined using an independently triggered sample of events. At lower en- 
ergies the triggers based on LAr information are supplemented by using additional information 
from the tracking detectors. In 1998 the LAr calorimeter electronics were upgraded in order to 
trigger scattered electrons with energies as low as 6 GeV, the minimum value considered in this 
analysis. This gives access to the high y kinematic region. For electron energies of 6 GeV, the 
overall trigger efficiency is 96% for the e+p data set and 90% for the earlier e'p data set. 

The characteristic feature of CC events is a large missing transverse momentum P^***, which 
is identified at the trigger level using the LAr calorimeter vector sum of "trigger towers", i.e. 
groups of trigger regions with a projective geometry pointing to the nominal interaction ver- 
tex. At low P™'*'^, the efficiency is enhanced by making use of an additional trigger requiring 
calorimeter energy in association with track information from the inner tracking chambers. For 
the minimum P^*** of 12 GeV considered in the analysis the efficiency is 60%, rising to 90% 
for P^''' of 25 GeV. In terms of the efficiency is 79% at 300 GeV^ and increases to 98% 
at 3 000 GcV^. These efficiencies are determined from the data using a sample of NC events 
in which all information from the scattered lepton is suppressed, the so-called pseudo-CC sam- 
ple. The trigger energy sums are then recalculated for the remaining hadronic final state. This 
sample also provides a useful high statistics cross check of further aspects of the CC analysis. 
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3.2 Simulation Programs 



Simulated DIS events are used in order to determine acceptance corrections. DIS processes are 
generated using the DJANGO [12] Monte Carlo (MC) simulation program, which is based on 
LEPTO [13] for the hard interaction and HERACLES [14] for single photon emission off the 
lepton line and virtual EW corrections. LEPTO combines 0{as) matrix elements with higher 
order QCD effects using the colour dipole model as implemented in ARL\DNE [15]. The 
JETSET program is used to simulate the hadronisation process [16]. In the event generation the 
DIS cross section is calculated with the PDFs of [17]. The simulated cross section is reweighted 
using a NLO QCD fit (HI 97 PDF fit) to previous data [1]. 

The detector response to events produced by the generation programs is simulated in detail 
using a program based on GEANT [18]. These simulated events are then subjected to the same 
reconstruction and analysis chain as the data. 

The dominant photoproduction background processes are simulated using the PYTHIA [19] 
generator with leading order PDFs for the proton and photon taken from [20]. Further back- 
ground from QED-Compton scattering, lepton pair production via two-photon interactions, 
prompt photon production and heavy gauge boson (W^, Z°) production are included in the 
background simulation. Further details are given in [1]. 

3.3 Polar Angle Measurement and Energy Calibration 

In the neutral current analysis the polar angle of the scattered electron (9^) is determined using 
the position of its energy deposit (cluster) in the LAr calorimeter, together with the interaction 
vertex reconstructed with tracks from charged particles in the event. The relative alignment of 
the calorimeter and tracking chambers is determined using a sample of events with a well mea- 
sured electron track, using information from both the CJC and the z drift chambers. Minimisa- 
tion of the spatial discrepancy between the electron track and cluster allows the LAr calorimeter 
and the inner tracking chambers to be aligned. The residual discrepancy then determines the 
systematic uncertainty on the measurement of 9e, which varies from 1 mrad for 9e > 135° to 
3mradfor^e < 120°. 

The calibration of the electromagnetic part of the LAr calorimeter is performed using the 
method described in [1]. Briefly, the redundancy of the detector information allows a predic- 
tion of the scattered electron energy (ED to be made based on the electron beam energy (Eg), 
the polar angle measurement of the scattered electron and the inclusive polar angle (7/j) [1] 
of the hadronic final state. This prediction of the double angle (DA) kinematic reconstruction 
method [21] is then compared with the measured electromagnetic energy, allowing local cali- 
bration factors to be determined in a finely segmented grid in z and 0. The calibration procedure 
is also performed on the simulated data. The final calibration is obtained by application of a 
further small correction determined from simulation, which accounts for small biases in the 
reconstruction of 7/^. The calibration is cross checked using independent data samples from 
QED-Compton scattering and two-photon e+e^ pair production processes. The total system- 
atic uncertainty on the absolute electromagnetic energy scale varies from 0.7% in the backward 
part of the calorimeter to 3% in the forward region, where statistics are limited. 
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Figure 1: Mean values of (a) Pr^h/ PT,e as a function of P^ g and (b) yh/VDA as a function of 
for neutral current data (solid points) and Monte Carlo (MC) simulation (open points) for 
7/i > 15° and 12 GeV < Pt,/i < 50 GeV. The curves correspond to a ±1% variation around the 
simulation. 



The hadronic final state is measured using energy deposits in the LAr and SPACAL calorime- 
ters supplemented by low momentum tracks. Isolated low energy calorimetric deposits are 
classified as noise and excluded from the analysis. The response of the detector to hadrons 
is calibrated by requiring transverse momentum conservation between the precisely calibrated 
scattered electron and the hadronic final state in NC events as described in [1]. The electron 
transverse momentum is defined as Px.e = \/v\,e + P^,e- The hadronic transverse momentum 
is determined from P^^ = a/ (Xli Vx,-i)^ + (Xli Vy,i)^^ where the summation is performed over 
all hadronic final state particles i, assuming particles of zero rest mass. 

Detailed studies and cross checks of the hadronic response of the calorimeter using the en- 
larged data sample have led to an improved understanding of the hadronic energy measure- 
ment. The calibration procedure is cross checked by requiring energy-momentum conserva- 
tion, E — Pz = (-Eg — Pz,e) + {Eh — Pz,h) = with Eh — Pz,h being the contribution of all 
hadronic final state particles Ylii^i — Pz,i)- In addition, the reference scale may be taken from 
the double angle method prediction rather than from the scattered electron. These studies have 
allowed the systematic uncertainty on the hadronic scale to be reduced with respect to previous 
measurements [1,3] in the region 12 GeV < Pt,/i < 50 GeV and 7/1 > 15°. The uncorrelated 
part (see section 3.6) of the hadronic scale uncertainty is reduced to 1% from 1.7% previously. 

Fig. 1 demonstrates the quality of the hadronic calibration in the stated region of Pt^h and 7/j, 
showing the level of agreement between data and simulation after the calibration procedure. 
In fig. 1(a) the mean value of the ratio PT,h/PT,e is shown. In fig. 1(b) the inelasticity yh, de- 
fined from the hadron reconstruction method [22] as yh = {Eh — Pz,h)/'^Ee, is compared with 
the DA variable yoA- In this analysis, it is the relative difference between data and simulation 
that is relevant and good agreement is found to within 1%. In addition a 1% correlated uncer- 
tainty is considered, accounting for possible remaining biases in the reference scale used for the 
calibrations. 
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3.4 Neutral Current Measurement Procedure 



Events from inelastic ep interactions are required to have a well defined interaction vertex to 
suppress beam-induced background. High Q"^ neutral current events are selected by requiring a 
compact and isolated energy deposit in the electromagnetic part of the LAr calorimeter^. The 
scattered electron is identified as the cluster of highest transverse momentum. In the central 
detector region, 9 > 35°, the cluster has to be associated with a track measured in the inner 
tracking chambers. 

As mentioned earlier, energy-momentum conservation requires E — — 2Ee. Restricting the 

measured E — P^to be greater than 35 GeV thus considerably reduces the radiative corrections 
due to initial state bremsstrahlung, where photons escape undetected in the backward direction. 
It also suppresses photoproduction background in which the scattered electron is lost in the 
backward beam-pipe and a hadron fakes the electron signal in the LAr calorimeter. Since the 
photoproduction background contribution increases with y, the analysis is separated into two 
distinct regions where different techniques are employed to suppress this background. The 
nominal analysis is restricted to < 0.63 for 90 GeV^ < Ql < 890 GcV^ and y^. < 0.90 for 
Ql > 890 GeVl This limits the minimum E'g to 11 GeV. The high-y analysis is performed 
for E'^> Q GeV, 0.63 < ye < 0.90 and 90 GeV^ < Ql < 890 GeVl Here ye and Ql are 
reconstructed using the scattered electron energy and angle, the so-called electron method. 

The NC kinematics in the nominal analysis are reconstructed using the eE method [23], which 
uses E^, 9e and Eh — Pz,h and has good resolution and small sensitivity to QED radiative 
corrections over the accessible phase space. In the high-y analysis the electron method gives 
the best resolution and is used to define the event kinematics. 

The nominal data sample consists of about 185 000 events. The comparison of the data and 

the simulation is shown in fig. 2 for the scattered electron energy and polar angle spectra and 
the distribution of E — P^, which are used in the reconstruction of x and Q^. All distributions 
are well described by the simulation, which is normalised to the luminosity of the data. In 
the nominal analysis the small photoproduction contribution is statistically subtracted using the 
background simulation. 

In the high-|/ analysis, the photoproduction background plays an increasingly important role, as 
low energies of the scattered electron are accessed. For this analysis, the calorimeter cluster of 
the scattered electron is linked to a well measured track having the same charge as the electron 
beam. This requirement removes a sizeable part of the background, where tt" 77 decays give 
rise to fake scattered electron candidates. The remaining background from photoproduction was 
estimated from the number of data events in which the detected lepton candidate has opposite 
charge to the beam lepton. This background is statistically subtracted assuming charge symme- 
try. The charge symmetry is determined to be 0.99 ± 0.07 by measuring the ratio of wrongly 
charged fake scattered lepton candidates in e~^p and e~p scattering, taking into account the dif- 
ference in luminosity. The charge symmetry is cross checked using a sample of data events in 
which the scattered electron is detected in the electron tagger and a systematic uncertainty of 
10% on the charge symmetry is assigned. Further details are given in [24, 25]. 

^Local detector regions are removed where the cluster of the scattered electron is not expected to be fuUy 
contained in the calorimeter, or where the trigger is not fuUy efficient. 
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Figure 2: Distributions of E'^ for (a) > 150 GeV^ and (b) > 5 000 GeV^ (c) 9^ and 
(d) E — for e^p data (solid points) and Monte Carlo (MC) simulation (open histograms) in 
the nominal analysis. The shaded histograms show the simulated background (bg) contribution, 
dominated by photoproduction. 



In total about 24 000 e^p events and 5 000 e^p events are selected in the high-y analysis. 
Figs. 3(a)-(c) show the scattered lepton energy spectrum, the polar angle distribution and the 
E — Pz spectrum for both the e^p and e~p data sets after background subtraction. The simula- 
tion, normalised to the luminosity of the data, provides a good description of these distributions. 
In fig. 3(d) the energy spectra of wrong charge lepton candidates in the data sets are shown. 
Good agreement is observed when the e^p data are normalised to the luminosity of the e~p 
data set. 



3.5 Charged Current Measurement Procedure 

The selection of charged current events requires a large missing transverse momentum, P^**^ = 
PT,h > 12 GeV, assumed to be carried by an undetected neutrino. In addition the event must 
have a well reconstructed vertex as for the NC selection. The kinematic variables and Q\ are 
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Figure 3: Distributions of (a) E'^, (b) 9^ and (c) E—P^ fore^p data (solid points), e^p data (open 
points) and Monte Carlo (MC) simulation (histograms) in the high-y analysis after background 
subtraction (see text). In (c) the E — cut is not applied, but is indicated by the dashed 
vertical line. Shown in (d) is the energy distribution of wrongly charged lepton candidates in 
background events. In (d) the e+p data have been normalised to the luminosity of the e'p data 
set. 



determined using the hadron kinematic reconstruction method [22]. In order to restrict the mea- 
surement to a region with good kinematic resolution the events are required to have yh < 0.85. 
In addition the measurement is confined to the region where the trigger efficiency is > 50% 
by demanding yh > 0.03. The ep background is dominantly due to photoproduction events in 
which the electron escapes undetected in the backward direction and missing transverse mo- 
mentum is reconstructed due to fluctuations in the detector response or undetected particles. 
This background is suppressed using the ratio Vap/Vp and the difference in azimuth between 

— * 

Pr,/i as measured in the main detector and the Plug calorimeter, A^/^ pi^g [3]. The quantities Vp 

— * 

and Vap are respectively the transverse energy flow parallel and anti-parallel to PT,h, the trans- 
verse momentum vector of the hadronic final state. The residual ep background is negligible for 
most of the measured kinematic domain, though it reaches 15% at the lowest and the highest 
y. The simulation is used to estimate this contribution, which is subtracted statistically from the 
CC data sample with a systematic uncertainty of 30% on the number of subtracted events. The 
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non-ep background is rejected as described in [1] by searching for event topologies typical of 
cosmic ray and beam-induced background. For further details see [26, 27]. 

After all selection criteria are applied, the final CC data sample contains about 1 500 events. The 
data and simulation are compared in fig. 4 for the Pt./i and Eh — Pz,h spectra, which are directly 
used in the reconstruction of the kinematic variables y and Q^. In both cases the simulation 
gives a reasonable description of the data. 

HI Charged Current Data 




20 40 60 80 100 10 20 30 40 



Px^h /GeV E,-P^, /GeV 

Figure 4: Distributions of (a) Pt^h ^.nd (b) — Pz,h for CC data (solid points) and Monte Carlo 
(MC) simulation (open histograms). The shaded histograms show the simulated background 
(bg) contribution, dominated by photoproduction. 



3.6 Cross Section Measurement and Systematic Uncertainties 

For both the NC and CC analyses the selected event samples are corrected for detector accep- 
tance and migrations using the simulation and are converted to bin-centred cross sections. The 
bins in the (x, Q"^) plane are defined as in refs. [1,3], based on the consideration of the detector 
resolution and event statistics. The bins used in the measurement are required to have values 
of stability and purity^ larger than 30%. This restricts the range of the NC measurements to 
y ^ 0.005. The QED radiative corrections (A^^^^p are defined in [1] and were calculated 
using the program HERACLES [14] as implemented in DJANGO [12]. 

The systematic uncertainties on the cross section measurements are presented in tables 6-13. 
They are split into bin-to-bin correlated and uncorrelated parts. All the correlated systematic 
errors are found to be symmetric to a good approximation and are assumed so in the following. 
The total systematic error is formed by adding the individual errors in quadrature. 

The correlated and uncorrelated systematic errors are discussed briefly below (see refs. [24-28] 
for more details). In addition, there is a global uncertainty of 1.5% and 1.8% on the luminosity 

^The stability (purity) was defined in [1] as the number of simulated events which originate from a bin and 
which are reconstructed in it, divided by the number of generated (reconstructed) events in that bin. 
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measurement for the e+p and e p data respectively, of which 0.5% is common to both (see 
section 4.1). 

• The total uncertainty on the electron energy scale is 1% if the z impact position of the 
electron at the calorimeter surface {zimp) is in the backward part of the detector {zimp < 
—150 cm), 0.7% in the region —150 cm < zi^p < 20 cm, 1.5% for 20 < Zi^p < 100 cm 
and 3% in the forward part (zi^p > 100 cm). The correlated part of the total uncertainty 
comes mainly from the possible bias of the calibration method and is estimated to be 0.5% 
throughout the LAr calorimeter. It results in a correlated systematic error on the NC cross 
section which is typically below 1%, increasing at low y to ~ 3% for Q^^l 000 GeV^ 
and ~ 8% for larger Q'^. 

• The correlated uncertainty on the electron polar angle is 1 mrad, 2 mrad and 3 mrad for 
9e > 135°, 120° < Oe < 135° and 9^ < 120°, respectively. This leads to a typical 
uncertainty on the NC reduced cross section of less than 1%, increasing up to ~ 5% at 
high X. 

• A0.5%(1%) uncorrected error originates from the electron identification efficiency in 

the NC nominal (high-y) analysis for Zimp < —5 cm. For Zimp > —5 cm the uncertainty is 
increased to 2%. The precision of this efficiency is estimated using an independent track 
based electron identification algorithm, limited for Zi^p > — 5 cm by statistics. 

• A 0.5% uncorrelated error is assigned on the efficiency of the scattered electron track- 
cluster link requirement in the NC nominal analysis. In the high-y analysis this is in- 
creased to 1%. 

• An uncorrelated 1% uncertainty on the hadronic energy measured in the LAr calorimeter 
is assigned for the region 12 GeV < Pt^h < 50 GeV and 7/j > 15°. Outside this region 
the uncertainty is increased to 1.7%. In addition, a 1% correlated component to the un- 
certainty is added in quadrature, originating from the calibration method and from the 
uncertainty on the reference scale. This yields a total uncertainty of 1.4% and 2% for the 
two regions respectively. 

• The uncertainty on the hadronic energy scale of the SPACAL calorimeter is 5%. The 
uncertainty on the hadronic energy measurement due to the inclusion of low momen- 
tum tracks is obtained by shifting their contribution by 3%. The influence on the cross 
section measurements from these sources is small compared with that from the corre- 
lated uncertainty from the LAr calorimeter energy scale. The three contributions (LAr, 
SPACAL, tracks) are thus combined, resulting in a single correlated hadronic error from 
the hadronic energy measurement, which is given in the tables. The corresponding error 
on the NC and CC cross sections is typically < 1%, but increases at low y to ~ 5%. 

• A 25% uncertainty is assigned on the amount of energy in the LAr calorimeter attributed 
to noise, which gives rise to a sizeable correlated systematic error at low y, reaching 
~ 10% at X = 0.65 and Q^<2 000 GcV^ in the NC measurements. 

• In the CC analysis the correlated uncertainties due to the cuts against photoproduction 
on Vap/Vp and A0/j.piug are only significant at high y, low Q'^ and low Pt,/i, reaching a 
maximum of ~ 7%. 
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• In the CC and the NC nominal analyses the photoproduction background is estimated 
from simulation. A 30% correlated uncertainty on the subtracted photoproduction back- 
ground is determined from a comparison of data and simulation for a phase space region 
dominated by photoproduction background. This results in a systematic error of typically 
<1%. 

• In the NC high-y analysis the photoproduction background is estimated directly from 
the data by using wrongly charged (fake) scattered lepton candidates, which leads to a 
10% correlated uncertainty on the subtracted photoproduction background. The resulting 
uncertainty on the measured cross sections is 1% or less. 

• A 0.3% uncorrelated error is considered on the trigger efficiency in the NC nominal anal- 
ysis and 2 — 6% in the CC analysis. For the NC high-y analysis the uncertainty on the 
cross section is ~ 2% at low Q^, decreasing to 0.6% at the highest in the analysis. 

• An uncorrelated error of 1% (NC) and 3% (CC) is estimated on the QED radiative cor- 
rections by comparing the radiative corrections predicted by the Monte Carlo program 
(DJANGO) with those calculated from HECTOR and EPRC [29]. The error also ac- 
counts for a small missing correction in DJANGO due to the exchange of two or more 
photons between the electron and the quark lines. 

• A 3% uncorrelated error is assigned on the event losses due to the non-ep background 
finders in the CC analysis, estimated from pseudo-CC data (see section 3.1). 

• A 2% uncorrelated error (5% for y < 0.1) on the vertex finding efficiency for CC events 
is estimated using pseudo-CC data. 

Overall a typical total systematic error of about 3% (6%) is reached for the NC (CC) double 
differential cross section. This precision has been achieved through detector improvements 
for triggering and a better understanding of the hadronic response of the detector, the electron 
identification and its angular measurement. 



4 QCD Analysis 

The cross section data presented here, together with the low precision data [30] and high 

e'^p data [1,3] previously published by the HI Collaboration, cover a huge range in 
and X. The improved accuracy now available allows the predictions of perturbative QCD to be 
tested over four orders of magnitude in Q"^ from about 1 GeV^ to above 10^ GeV^, and x from 
below 10~^ to 0.65. The measurements of NC and CC e'^p scattering cross sections provide 
complementary sensitivity to different quark distributions and the gluon distribution, xg{x, Q^). 
This is used to determine the sum of up-type xU, of down-type xD and of their anti-quark-type 
xU and xD distributions, employing only HI inclusive cross section data. 

With the current beam energies, the HERA collider data do not give access to the large x region 
of deep inelastic scattering at medium (~ 100 GeV^). Complementary information on quark 



15 



distributions in this region is provided by fixed-target lepton-proton data. Lepton-deuteron scat- 
tering data, which provide further constraints on the PDFs, are not yet available from HERA. 
Therefore, in the subsequent analysis, the HI data are also combined with the precise BCDMS 
muon-proton and muon-deuteron scattering data and the results are compared for cross checks 
of the PDFs obtained from the analysis of the HI data alone. 

4.1 Ansatz 

Traditionally, QCD analyses of inclusive deep inelastic scattering cross section data have used 
parameterisations of the valence quark distributions and of a sea quark distribution, imposing 
additional assumptions on the flavour decomposition of the sea [1, 31-33]. The neutral and 
charged current cross section data presented here, however, are sensitive to four combinations 
of up- and down-type (anti-)quark distributions which, for less than the bottom quark pro- 
duction threshold, are given above in eq. 10. Working in terms of these combinations weakens 
the influence of necessary assumptions on the flavour decomposition of the sea in the fit. The 
valence quark distributions are obtained from 

xu^ ^x{U-U) , xdy = X {D -D) (12) 

and are not fitted directly. 

In the QPM, the charged current structure function terms 0^,^ are superpositions of the distri- 
butions given in eq. 10 according to 

^xU + {l- yfxD, 4>-cc = xC/ + (1 - yfxD . (13) 

The neutral current structure function terms (p^Q are dominated by the electromagnetic structure 
function F2 , which can be written as 

F2^^x{U + U) + '^-x{D + D) . (14) 

In the high neutral current data, complementary sensitivity is obtained from the interference 
structure function xF^^ = x[2{U — U) + {D — D)]/3, but still higher luminosity is required to 
exploit this for a dedicated determination of the valence quarks. 

In the fit to the HI and BCDMS data, the isoscalar nucleon structure function is determined 
by the singlet combination of parton distributions and a small contribution from the difference 
of strange and charm quark distributions, 

Fi" ^^x{U + U + D + D) + lx{c + c-s-s) . (15) 

lo D 

The nucleon data obtained from the BCDMS muon-deuteron cross sections are measured for 
X > 0.07. For these data nuclear corrections are applied following [34]. In eq. 15 the charm 

and strange quark distributions occur explicitly and may be constrained using experimental data 
as provided by HI and ZEUS on the charm contribution to F2 [35,36] and from NuTeV on the 
strangeness content of the nucleon [37]. The analysis of the HI data, however, is rather in- 
sensitive to these quark distributions. They are assumed to be fixed fractions of the up- and 
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down-type quark distributions respectively at the initial scale of the QCD evolution (see sec- 
tion 4.2). 



The analysis is performed in the MS renormalisation scheme using the DGLAP evolution equa- 
tions [38] at NLO [39] . The structure function formulae given here are thus replaced by integral 
convolutions of coefficient functions and PDFs. An approach is used whereby all quarks are 
taken to be massless, including the charm and bottom quarks, which provides an adequate de- 
scription of the parton distributions in the high kinematic range of the new data presented 
here. The bottom quark distribution, xh, is assumed to be zero for < ml where m;, is the 
bottom quark mass. 

Fits are performed to the measured cross sections calculating the longitudinal structure func- 
tions to order a1 and assuming the strong coupling constant to be equal to as{M'^) = 0.1185 [7]. 
All terms in eqs. 1 and 7 are calculated, including the weak corrections, A^'^^^. The analysis 
uses an x space program developed inside the HI collaboration [40], with cross checks per- 
formed using an independent program [41]. In the fit procedure, a function is minimised 
which is defined in [30]. The minimisation takes into account correlations of data points caused 
by systematic uncertainties allowing the error parameters (see table 2), including the relative 
normalisation of the various data sets, to be determined by the fit. The fit to only HI inclusive 
cross section data, termed HI PDF 2000, uses the data sets as specified in table 1. The table 
additionally lists the BCDMS data used in a further fit for comparison with the HI PDF 2000 
fit. 



data set 


process 


X range 




range 


S'- 


ref. 


comment 










(GeV^) 


(GeV2) 


(%) 






HI minimum bias 97 


e+p NC 


0.00008 


0.02 


1.5 


12 


1.7 


[30] 


Vi = 301 GeV 


HI low Q2 96 _ 97 


e+p NC 


0.000161 


0.20 


12 


150 


1.7 


[30] 


= 301 GeV 


HI high Q2 94 _ 97 


e+pNC 


0.0032 


0.65 


150 


30 000 


1.5 


[1] 


V5 = 301 GeV 


HI high Q2 94 _ 97 


e+pCC 


0.013 


0.40 


300 


15 000 


1.5 


[1] 


= 301 GeV 


HI high (52 98 - 99 


e-pNC 


0.0032 


0.65 


150 


30 000 


1.8 


[3] 


Vs = 319 GeV 


HI high (32 98 _ 99 


e-pCC 


0.013 


0.40 


300 


15 000 


1.8 


[3] 


Vi = 319 GeV 


HI high Q2 98 _ 99 


e-pNC 


0.00131 


0.0105 


100 


800 


1.8 


this rep. 


Vi = 319 GeV; high-y data 


HI high Q2 99 _ 00 


e+pNC 


0.0032 


0.65 


150 


30 000 


1.5 


this rep. 


= 319 GeV; incl. high-j/ data 


HI high Q2 99 _ 00 


e+pCC 


0.013 


0.40 


300 


15 000 


1.5 


this rep. 


Vi = 319GcV 


BCDMS-p 


HpNC 


0.07 


0.75 


7.5 


230 


3.0 


[4] 


require yp, > 0.3 


BCDMS-D 


/iDNC 


0.07 


0.75 


7.5 


230 


3.0 


[4] 


require y/j, > 0.3 



Table 1: Data sets from HI used in the HI PDF 2000 fit and from BCDMS /i-proton and fx- 
deuteron scattering used in the Hl+BCDMS fit. As for the previous HI QCD analysis [30], the 
original BCDMS data are used at four different beam energies imposing the constraint > 0.3. 
The inelasticity was defined using BCDMS beam energies. The normalisation uncertainties 
of each data set (S^) are given as weJJ as the kinematic ranges in x and Q^. The uncertainty 
5^ includes a common error of 0.5% for the HI data sets (see text). The nominal analysis and 
high-y analysis do not overlap in kinematic coverage (see section 3.4). 

The correlated systematic uncertainties for the HI cross section measurements may be corre- 
lated across data sets as well as between data points, since they may arise from the same source. 
They are thus not treated independently in the QCD analysis presented here. The relationship 
between the error sources as used in the fitting procedure is summarised in table 2 for each of 
the eight correlated systematics considered. This leads to 18 independent error parameters. In 
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addition, all HI quoted luminosity uncertainties have a common contribution of 0.5% arising 
from the theoretical uncertainty on the Bethe-Heitler cross section. This common contribution 
has been taken into account in the QCD analysis. 



data set 




process 


()^ 


S^' 




()'' 




()^' 




- Q 

()^ 


HI minimum bias 97 


e+p NC 


CI 


El 


91 


hi 


Nl 


Bl 






HI low Q2 96 - 


- 97 


e+p NC 


C2 


El 


61 


hi 


Nl 


Bl 






HI high g2 94 


-97 


e+pNC 


C3 


E2 


92 


h2 


Nl 


B2 






HI high g2 94 


-97 


e+pCC 


£3 






h2 


Nl 


B2 


VI 




HI high 98 


-99 


e~pNC 


£4 


E2 


93 


h2 


Nl 


B2 




SI 


HI high g2 98 


-99 


e-p CC 


£4 






h2 


Nl 


B2 


V2 




HI high g2 99 


-00 


e+p NC 


£5 


E2 


93 


h2 


Nl 


B2 




SI 


HI high g2 99 


-00 


e+pCC 


C5 






h2 


Nl 


B2 


V2 





Table 2: Treatment of the correlated systematic error sources for the HI data sets used in the 
fits. For each of the eight correlated systematic error sources, one or more parameters are 
included in the QCD fit procedure. The sources considered are due to the luminosity uncertainty 
(S^), the electron energy uncertainty (5^), the electron polar angle measurement uncertainty 
(5^), the hadronic energy uncertainty (5^), the uncertainty due to noise subtraction (6^), the 
photoproduction dominated background simulation error (6^), the uncertainty due to the cuts 
against photoproduction in the CC analysis (6^) and the error on the charge symmetry in the 
high-y analysis (S^). The table entries indicate the correlation of the error sources across the 
HI data sets. For example, the uncertainty due to the noise subtraction is the same for all data 
sets leading to one common parameter in the fit (Nl), whereas the electron energy uncertainty 
has two independently varying parameters (El and E2 ) for the HI JVC data sets only. 



4.2 Parameterisations 

The initial parton distributions, xP — xg, xU, xD, xU, xD, are parameterised at g^ = Ql in 
the following general form 

xP{x) = Apx^^{l - x f''[l + Dpx + Epx'^ + Fpx-^ + Gpx% (16) 

The QCD analysis requires choices to be made for the initial scale (Ql) and the minimum g^ 
of the data considered in the analysis (g^m)- Variations of both Ql and g^jn are studied. As 
in [30] Ql is chosen to be 4 GeV^ and g^j„ = 3.5 GeV^. Reasonable variations of these choices 
are considered as part of the model uncertainties on the parton distributions (section 5.3). 

The general ansatz, eq. 16, represents an over-parameterisation of the data considered here. The 
specific choice of these parameterisations is obtained from saturation of the x^'- additional 
parameter D, E,ForG is considered only when its introduction significantly improves the x^- 
The appropriate number of parameters also depends on the data sets included in the fit. The 
HI data requires less parameters than the combined HI and BCDMS data due to the precise 
BCDMS proton and deuteron data in the large x region, where the cross section variations with 
X are particularly strong. 
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The fit to the HI data uses the following functional forms 

xg{x) = AgX^'{l - x)'^' ■ [1 + Dgx] 

xU{x) = Aux^^{l - xf^ ■ [1 + Dux + Fux^] 

xD{x) = Adx^''{1 - ■ [1 + Ddx] (17) 

xU(x) = Ajjx^{l - x)^ 

xD{x) = A-^'^il- x)'^T^ , 

in which the number of free parameters are further reduced using the constraints and assump- 
tions detailed below. 

The number of parameters required by the fit for the different parton distributions follows the 
expectation. A high x term Egx'^ is not needed in the gluon distribution, since at large x the 
scaling violations are due to gluon bremsstrahlung, i.e. are independent of the gluon distribution. 
The xU and xD distributions require more parameters than the anti-quark distributions xU 
and xD because the former are a superposition of valence and sea quarks, in contrast to the 
latter. Due to the different electric charges, = 4e^, and the y dependence of the charged 
current cross section, the data are much more sensitive to the up quark than to the down quark 
distributions. Thus less parameters are needed for xD than for xU. 

A number of relations between parameters can be introduced naturally in this ansatz. At low 
X the valence quark distributions are expected to vanish and the sea quark and the anti-quark 
distributions can be assumed to be equal. Thus the low x parameters Ag and Bg are required to 
be the same for xU, xU and for xD, xD. In the absence of deuteron data from HERA there is 
no distinction possible of the rise towards low x between xU and xD. Thus the corresponding 
B parameters are required to be equal, i.e. Bu — Bd — Bjj — B^ = Bg. Further constraints 
are the conventional momentum sum rule and the valence quark counting rules. 

The ansatz presented above allows the quark distributions xU, xD, xU, xD to be determined. 
Further disentangling the individual quark flavour contributions to the sea is possible only with 
additional experimental information and/or assumptions. Assuming that the strange and charm 
sea quark distributions xs and xc can be expressed as ^-independent fractions fs and fc of xD 
and xU at the starting scale ofQl = A GeV^ (see table 5), a further constraint is used in the fit: 
Ajj — A-^ ■ (1 — fs) /(I — fc), which imposes that d/u — > 1 as x — > 0. 

The total number of free parameters of the five parton distributions is thus equal to 10 in the fit 
to the HI data. The value is hardly improved by including any half integer power of x. The 

parametric form of eq. 17 is also found starting from an alternative polynomial in x^, which 
includes half integer powers up to x^/^ . The addition of the large x BCDMS jip and jjiD data 
leads to two additional terms, Gux^ and FdX"^, in the polynomials. 



5 Results 

5.1 NC and CC Cross Sections dcr /dQ^, dcr /dx and erg'*. 

The e^p single differential neutral current cross section da/dQ^ measured for y < 0.9 is 
shown in fig. 5(a). The data are compared with previous HI e+p measurements made at 
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= 301 GeV. The new cross sections are found to be higher than the measurement from 

94 — 97 as expected due to the increase in centre-of-mass energy. Both cross sections, falling by 
over six orders of magnitude for the measured region between 200 GcV^ and 30 000 GcV^, 
are well described by the HI PDF 2000 fit. The error band represents the total uncertainty as 
derived from the QCD analysis by adding in quadrature the experimental and model uncer- 
tainty. The experimental uncertainty on the predicted cross sections is significantly larger than 
the model uncertainty, which is discussed in section 5.3. Fig. 5(b) shows the ratios of the mea- 
surements to the corresponding Standard Model expectation determined from the HI PDF 2000 
fit. Note that in this lower figure the HI data are scaled by the normalisation shift imposed by 
the QCD fit given in table 3. The new data are given in table 6. 

The dependence of the charged current cross section from the 99 — 00 data is shown in 
fig. 6(a). For consistency with the NC cross sections, the data are presented in the range |/ < 0.9, 
after correction"^ for the kinematic cuts 0.03 < y < 0.85 and Pt^h > 12 GeV (section 3.5). 
The data are compared with the previous measurement taken at lower centre-of-mass energy. 
The ratios of data to expectations are shown in fig. 6(b) together with the Standard Model 
uncertainty. Again in this lower figure the HI data are scaled by the normalisation shift imposed 
by the QCD fit, given in table 3. The two data sets agree well with each other, though the new 
data have a tendency to be higher than the fit result at high Q^. The data are listed in table 7. 

Fig. 7 shows the Q'^ dependences of the NC and CC cross sections representing the total e+p 
and e~p data sets taken at HERA-I. The e+p data have been combined after scaling the 94 — 97 
data to = 319 GeV, using the HI PDF 2000 fit and the procedure described in [27]. At low 
the NC cross section exceeds the CC cross section by more than two orders of magnitude. 
The sharp increase of the NC cross section with decreasing Q'^ is due to the dominating photon 
exchange cross section with the propagator term oc l/Q^. In contrast the CC cross section 
(~ [M^/((5^ + M^)]^) approaches a constant at low Q^. The CC and NC cross sections are of 
comparable size at ~ 10^ GeV^, where the photon and Z° exchange contributions to the NC 
cross sections are of similar size to those of exchange. These measurements thus illustrate 
unification of the electromagnetic and the weak interactions in deep inelastic scattering. While 
the difference in NC cross sections between e^p and e^p scattering is due to interference, 
the difference of CC cross sections arises from the difference between the up and down quark 
distributions and the less favourable helicity factor in the e+p cross section (see eq. 13). 

The single differential cross sections da/dx are measured for Q'^ > 1 000 GeV^ for both NC 
and CC and also for > 10 000 GeV^ in the NC case. The NC data are compared in fig. 8 with 

the previous HI e~^p measurement at ^/s = 301 GeV and the corresponding expectations from 
the fit. A similar comparison for the CC data is shown in fig. 9. Increases with ^/s are observed 
in both the NC and the CC cross sections, in agreement with the expectations. The fall in the 
cross sections at low x is due to the restriction y < 0.9. The measurements are summarised in 
tables 8-10. 

The total CC cross section has been measured in the region Q'^ > 1 000 GeV^ and y < 0.9 after 
applying a small correction factor of 1.03 for the y and Pr,/i cuts, determined from the HI PDF 
2000 fit. The result is 

(T^°^(e+p; ^ = 319 GeV) = 19.19 ± 0.61(stat.) ± 0.82(syst.) pb, 
"^The correction factors are given in table 7. 
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where the 1.5% normalisation uncertainty is included in the systematic error. This is to be 
compared with the value from the HI PDF 2000 fit (T^^(e+p) = 16.76±0.32 pb. The difference 
between the measurement and the fit is 2.3 standard deviations assuming the correlation of 
uncertainties between measurement and fit is negligible. An unbiased theoretical expectation 
for (T^cci^^P) i^^y be obtained by repeating the HI PDF 2000 fit but excluding the new 99 — 00 
CC data, which leads to 16.66 ± 0.54 pb. 

Additionally, the analysis has been performed on the 94 — 97 data set at the lower centre-of-mass 
energy, yielding 

y/s = 301 GeV) = 16.41 ± 0.80(stat.) ± 0.90(syst.) pb . 

This is to be compared with the cross section obtained from the HI PDF 2000 fit ^^^(e+p) = 
14.76 ± 0.30 pb. Assuming that the systematic uncertainties are fully correlated and part of the 
luminosity uncertainties are common (section 4.1), the 94 — 97 and 99 — 00 measurements are 
combined [27] yielding a value of 

(^cc(e^P; = 319 GeV) = 18.99 ± 0.52(stat.) ± 0.81(syst.) pb . 

5.2 NC and CC Double Differential Cross Sections 

The double differential NC reduced cross section, unc (defined in eq. 11), is shown in fig. 10 
for both the nominal and high-|/ 99 — 00 e+p data. In addition the new high-y 98 — 99 e~p data 
are presented. The data agree well with the expectations of the HI PDF 2000 fit, which are also 
shown^. The rise of the Standard Model DIS cross section towards low x (high y) departs from 
the monotonic behaviour of F2 due to the contribution of the longitudinal structure function Fl. 
This allows Fl to be determined in the high y region (section 5.5). 

In fig. 11 the e+p NC large x cross section data at ^/s — 319 GeV are compared with the data 
obtained previously [1] at y/s = 301 GeV. The two data sets are found to be in agreement 
with each other and with the HI PDF 2000 fit. Fig. 1 1 also shows the data from the recent HI 
measurement at lower Q"^ [30] and the fixed-target data from BCDMS [4] . The fit description 
of the BCDMS data, which are not used in the fit, is remarkably good except at very large 
X — 0.65. A similar observation has already been reported in [1, 30]. At the highest a 
decrease of the cross section is expected due to the negative 7Z interference in e+p scattering. 

In fig. 12 the reduced CC cross section, dec (defined in eq. 11), is shown for the new data and 
the data taken at lower energy between 1994 and 1997. These data are found to be compatible 
with each other considering the weak energy dependence of the reduced CC cross section. An 
extension of the x range for = 3 000 GeV^ and 5 000 GcV^ is achieved due to the improved 
trigger efficiency. The combined 94 — 00 result is compared in fig. 13 with the expectation from 
the HI PDF 2000 fit. Also shown is the expected contribution of the xD distribution, which 
dominates the cross section at large x. The HERA e+p CC data can thus be used to constrain 
the d quark distribution in the valence region. 

All double differential measurements together with the contributions of each of the major sys- 
tematic uncertainties are listed in tables 11-13. 

^The normalisation factors as determined by the QCD fit (table 3) are not applied to the data shown in the 
figure. 
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5.3 Fit Results 



In this section, the results of the QCD analysis are presented. The value for each data set 
is given in table 3 as well as the optimised relative normalisation as determined from the fit. 
The total value^ per degree of freedom (x^/ndf) is 540/(621 - 10) = 0.88. The NLO QCD 
fit requires the lowest data (HI minimum bias 97, < 12 GeV^) to be raised by 3.7%, 
corresponding to 2.3 standard deviations in terms of the normalisation uncertainty excluding the 
common error of 0.5% (see section 4.1). In contrast all data for Q'^ > 100 GeV^ are lowered, 
by at most 1.9%. It can not yet be decided whether this behaviour is due to inadequacies in the 
theory (e.g. the missing higher order terms in In Q"^) or experimental effects. 



data set 


process 


data points 


(unc. err.) 




(corr. err.) 


normalisation 


HI minimum bias 97 


e+p NC 


45 


37.5 




5.9 


1.037 


HI low g2 96 _ 97 


e+pNC 


80 


71.2 




1.3 


1.008 


HI high g2 94 _ 97 


e+pNC 


130 


89.7 




2.1 


0.981 


HI high g2 94 _ 97 


e+pCC 


25 


18.0 




0.4 


0.981 


HI high g2 98 - 99 


e-p NC 


139 


114.7 




1.0 


0.991 


HI high 98 - 99 


e-p CC 


27 


19.5 




0.7 


0.991 


HI high g2 99 _ 00 


e+pNC 


147 


142.6 




2.6 


0.985 


HI high Q2 99 _ 00 


e+pCC 


28 


32.4 




0.9 


0.985 


Total 


621 


540 





Table 3: For each data set used in the HI PDF 2000 Gt, the number of data points is shown, 
along with the contribution determined using the uncorrelated errors (unc. err). Each of the 
correlated error sources (see table 2) leads to an additional contribution [30], which is listed as 
(corr err). Also shown is the optimised normalisation of the data set as determined by the 
m. The HI NC 98 - 99 e'p and HI NC 99 - 00 e+p data include the high-y analyses. 

The parameters of the initial parton distributions are given in table 4 (see also [42]) and the 
distributions are shown in fig. 14. The inner error band describes the experimental uncertainty, 
while the outer band represents the experimental and model uncertainties added in quadrature. 

The experimental accuracy of the initial distributions is typically a few percent in the bulk of 
the phase space of the HI data. This accuracy has negligible dependence on but a strong 
dependence on x. The best precision is achieved for the xU quark distribution, which amounts 
to about 1% for x = 0.01 and reaches 3% and 7% at x = 0.4 and 0.65, respectively. The xD 
quark distribution is only determined with moderate precision as it is predominantly constrained 
by the CC e+p cross sections, which are still subject to limited precision. The corresponding 
uncertainties on xD at the three quoted x values are respectively ~ 2%, ~ 10% and ~ 30%. 

These uncertainties reflect the kinematic dependence and size of the measurement errors. How- 
ever the error size also depends significantly on the fit assumptions. If, for example, the con- 
straint between Ajj and A-^ on the low x behaviour of the anti-quark distributions is relaxed, 
the small uncertainty at low x — 0.01 is much increased to ~ 6% and ~ 20% respectively for 

*In the calculation of the x^, the assumption is made that the uncorrelated errors among different data points 
within one data set stay uncorrelated with the corresponding data points from an independent data set. 
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xU and xD. The measurement of the low x behaviour of up and down quarks and their possible 
distinction requires electron-deuteron data to be taken at HERA. 

The model parameter uncertainties on the PDFs are determined in a similar manner to [30] and 
the sources of uncertainty are specified in table 5. The model uncertainties are relatively small 
with respect to those from experimental sources except at small x and low Q^, where they reach 
~ 2% and ~ 6% respectively for xU and xD at x = 0.01 and = 10 GeVl 

Within the functional form considered (see eq. 16), the parameterisation given in eq. 17 is found 
uniquely. Possible variations within the Ax^ — 1 region of the parameter space do not lead to 
noticeably different distributions. Thus in this analysis no account is made of uncertainties 
due to the choice of parameterisations. A completely different ansatz, however, may well lead 
to different initial distributions, as seen, for example, in the complicated shape of xg chosen 
in [33]. The gluon distribution determined in this analysis is consistent with the distribution 
obtained previously by HI [30] if the effects of the different heavy flavour treatments are taken 
into account. 



p 


A 


B 


C 


D 


F 


xg 


0.0183 


-0.872 


8.97 


3450. 




xU 


0.112 


-0.227 


5.08 


48.0 


373. 


xD 


0.142 


-0.227 


4.93 


23.5 




xU 


0.112 


-0.227 


7.28 






xD 


0.142 


-0.227 


4.36 







Table 4: Parameters of the HI PDF 2000 fit to the HI data alone for the initial distributions 
at Ql = A GeV^. Equal parameter values reflect the constraints imposed by the fit (see sec- 
tion 4.2). The uncertainties and their correlations are available in [42]. 



source of uncertainty 


central value 


variation 


QLn (GeV^) 


3.5 


2.0 - 5.0 


Ql (GeV2) 


4.0 


2.0 - 6.0 


«s(M|) 


0.1185 


0.1165 - 0.1205 


fs, strange fraction of xD 


0.33 


0.25 - 0.40 


fc, charm fraction of xU 


0.15 


0.10 - 0.20 


rric (GeV) 


1.4 


1.2 - 1.6 


rub (GeV) 


4.5 


4.0 - 5.0 



Table 5: Model uncertainties considered in the QCD analysis. 



The full curve in fig. 14 is the result [43] of the fit to HI and BCDMS data, which gives a 
X^/ndf= 883/(1014 - 12) = 0.88. Excellent agreement of the PDFs between the two fits is 
observed. For large x, the high data of HI allow distinction between up and down flavours 
yielding results compatible with those from BCDMS proton and deuteron data. At low x only 
HERA data are available and thus the two fits are forced to be in agreement, apart from possible 
small effects due to sum rules. 
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The PDFs from the HI PDF 2000 fit are further compared in fig. 15 with recent results from 
the MRST [32] and CTEQ [33] groups for two values of at 10 GeV^ and 1 000 GcV^. The 
HI PDF 2000 fit is in remarkable agreement with the MRST and in particular the CTEQ anal- 
yses, given the many differences in terms of the data sets used, the assumptions made and the 
treatment of heavy flavours. 



5.4 Extraction of the Proton Structure Function F2 

The NC structure function term 07vc is obtained from the measured NC double differential cross 
section according to eq. 1. It is dominated by the structure function F2 in most of the kinematic 
range. The structure function F2 is extracted using 

= ^(1 + A^, + + A^J-^ . (18) 

Here the correction terms ^p,^ and A^?,^ account for the effects of exchange on F2 and xF-^ 
(eqs. 2-4) and originates from the longitudinal structure function F^. These corrections, 
shown in table 11, are determined using the HI PDF 2000 fit (see section 4). At high y and 
< 1 500 GeV^, A^^ is sizeable and the extraction of F2 in this Q"^ region is thus restricted 
to the kinematic range y < 0.6. It is extended to higher y at larger (> 2 000 GeV^) where 
the predicted contribution of Fl is small. 

The extracted F2 using the 99 — 00 data is presented in table 11. Fig. 16 shows the F2 data 
using the combined 94 — 97 and 99 — 00 high Q"^ e+j; data sets. Also shown in the figure are the 
recent HI F2 data at lower [30] and structure function data from BCDMS [4] and NMC [44]. 
The full HI data cover a range of four orders of magnitude in x and Q"^. The HI PDF 2000 fit 
provides a good description of the data over the whole region except for the BCDMS data at 
X = 0.65, as seen in fig. 11. The fit also gives a good description of the scaling violations 
observed in the measurements. 



5.5 Determination of the Longitudinal Structure Function Fl 

The structure function term 0jvc is used to determine Fl at y > 0.63 and < 890 GeV^. For 
statistical reasons, the measured cross sections in two neighbouring bins are combined, as- 
suming that the systematic uncertainties are fully correlated. The longitudinal structure function 
is then determined using the formula 

Y+F2TY-xF-i-<l>%c] , (19) 

for e^p scattering which, neglecting the small electroweak contributions in the region of this 
extraction, reduces to the expression 

= 1 [Y^F2 - <l>%c] ■ (20) 
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The extraction of relies upon the extrapolation of the fit result for into the high y region, 
that is, to larger for given x. In order to avoid a possible influence of the high y data region 
on this calculation, a dedicated NLO QCD fit (HI Low y fit) is performed to HI data with 
y < 0.35 only and the results are extrapolated using the DGLAP evolution equations. This 
method was introduced in [45]. 

Apart from the y range restriction, the HI Low y fit follows the same procedure as described 
in section 4. It results in a x^/ndf= 417/ (455 — 10) — 0.94 and agrees very well with the HI 
PDF 2000 fit over the full y range. The normalisation shifts of the data sets used are found to 
be within 1% of those from the HI PDF 2000 fit. 

In the extraction of the longitudinal structure function, the experimental cross sections are 
slightly modified using the results of the HI Low y fit for the renormalisation and small shifts 
from the correlated uncertainties common to the low y and the high y region. The combined 
HERA-I measurements of the structure function term 0jvc and the extracted values of F^ are 
listed in table 14. The statistical precision is due directly to the cross section measurement 
at high y. The systematic uncertainties arise from the measurement errors at high y and the 
model uncertainties related to the extrapolation of F2 from the low y to the high y region. The 
correlations in the systematic uncertainties between low and high y are taken into account. 

In fig. 17 the determinations of F^ at high Q"^ are shown for both the e^p and the e~p data sets. 
The results from both data sets are mutually consistent and are in agreement with the HI Low y 
fit prediction for Fl, based on the gluon distribution derived from the scaling violations of F2 
at larger x. The extreme values allowed for Fl {Fl = and Fl = F2) are clearly excluded by 
the data. A model independent measurement of Fl and the x dependence can, however, only 
be achieved with reduced beam energies at HERA. 

5.6 Measurement of the Generalised Structure Function xF^ 

At high Q'^, the NC cross section in e+p scattering has been observed to be significantly smaller 

than that in e~p scattering [3], confirming the Standard Model expectation of 7Z interference. 
These earlier HI data were used to obtain a first measurement of the generalised structure func- 
tion xF^ in the kinematic range 0.02 <x< 0.65 and 1 500 GeV^ < < 12 000 GeV^ [3]. A 
similar measurement has been reported recently by ZEUS [46]. 

Profiting from the enlarged statistics and the reduced systematic uncertainties, the previous 
measurement of xF^ [3] is updated here by using the same published e~p and the full e+p data 
obtained by HI at HERA-I. Fig. 18(a) shows the comparison of the e~'p and e'^p data for three 
different values at 1 500 GcV^ 5 000 GeV^ and 12 000 GeV^ together with the expectations 
determined from the HI PDF 2000 fit. The generalised structure function xF^, given in table 15, 
is obtained from 

^Fs = ^ - 'Pnc] (21) 

and is compared in fig. 18(b) with the calculation. Since at high x and low Q"^ the expected 
sensitivity to xF-^ is smaller than the luminosity uncertainty, the measurement is not performed 
in this region. The dominant contribution to xF^ arises from interference, which allows 
xF^^ to be extracted (eq. 4) according to xF^^ ~ —xF^iQ"^ + M|)/(aefi;(5^) by neglecting 
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the pure Z exchange contribution, which is suppressed by the small vector coupling v^. This 

structure function is non-singlet and has little dependence on Q^. This is illustrated in fig. 18(c). 
The measured xF]^"^ at these values can thus be averaged taking into account the small 
dependence. The two lowest x bins at a; = 0.020 and 0.032 are averaged as well. The 
averaged xF^^, determined for a value of 1 500 GeV^, is shown in table 15 and fig. 18(d) in 
comparison with the QCD fit result. The structure function xF^^ determines both the shape and 
magnitude of the valence quark distributions independent of the sea quark distributions. The 
calculation from the QCD fit, in which the parton densities in the valence region are principally 
constrained by the NC and CC cross sections rather than the difference between the NC 
cross sections, gives a good description of the measurement. The averaged structure function is 
integrated [47] over the measured x range, yielding 

/■0.650 

/ F^^{x, = 1 500GeV^)da; = 1.28 ± 0.17(stat.) ± 0.11(syst.) , 

Jo. 026 

which is in agreement with 1.06 ± 0.02, as predicted from the HI PDF 2000 fit. 

5.7 The Quark Distributions xu and xd at Large x 

The flavour composition of the proton at high x may be disentangled by exploiting the NC 
and CC cross section measurements. The e+p CC cross section at large x is dominated by the 
d quark contribution as may be inferred from fig. 13. Similarly the u distribution dominates 
the e~p CC and e^p NC cross sections at large x. Using data points for which the xu or xd 
contribution provides at least 70% of the cross section, as determined from the HI PDF 2000 
fit, the up and down quark distributions are determined locally, using the method described 
in [1,27]. The extraction relies on weighting the differential cross section measurement with the 
calculated local flavour contribution and is illustrated in fig. 19, where xu is the combined result 
from three independent extractions from the NC e'^p and CC e~p data and xd is determined 
from the CC e+p data only. This method is complementary to performing a QCD fit, since it is 
based on the local cross section measurements and is less sensitive to the parameterisations and 
dynamical assumptions used in the fit. 

The extracted xu and xd distributions are further compared in fig. 19 with the results of the 
HI PDF 2000 fit by subtracting xc and xs from the fitted xU and xD. The two determinations 
are in good agreement. They also compare well with the recent parameterisations from the 
MRST [32] and CTEQ [33] groups except for xu at large x = 0.65, where the results of MRST 
and CTEQ, being constrained mainly by the BCDMS data, yield a larger up quark distribution. 

6 Summary 

New measurements are presented of inclusive deep inelastic neutral and charged current scatter- 
ing cross sections at high momentum transfers > 100 GcV^ from recent e^p data recorded 
in 1999 and 2000 by the HI experiment at HERA. This analysis, together with previous analy- 
ses of the 1994 — 1997 e+p and 1998 — 1999 e^p data, completes the HI measurements of the 
inclusive cross sections at high Q"^ from the first phase of HERA operation. 
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The accuracy of the neutral current (NC) measurements presented here has reached the level of 
a few percent in the medium range of < 3 000 GcV^. The very high NC and charged 
current (CC) data are still limited by their statistical precision, which is expected to improve in 
the high luminosity phase of HERA. 

For both e+p and e^p scattering data, the region of very large inelasticity is explored, which 
allows a determination of the longitudinal structure function Ft(a;, Q^) for the first time in the 
large momentum transfer range, 110 GeV^ < < 700 GeV^. The observed interference of 
the photon and Z exchange, differing between e+p and e~p NC scattering at high Q^, is used 
to measure the structure function superseding the earlier measurement. 

The NC and CC cross sections in e^ip scattering are sensitive to the sums of up- and down-type 
quark and anti-quark distributions, xU , xD, xU and xD. Based on these quark distribution 
decompositions, a novel NLO QCD analysis is performed, resulting in a first determination of 
the partonic nucleon structure from inclusive DIS scattering data from HI alone. The experi- 
mental precision achieved in this analysis is about 3% and 10% respectively for xU and xD at 
X — OA. The resulting parton distributions are found to be in agreement with those obtained in 
an analysis also including the BCDMS muon-nucleon data at large x. The QCD analysis leads 
to a good description of all the fitted NC and CC cross section data and of the derived structure 
functions over more than four orders of magnitude in x and Q'^. 
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Figure 5: (a) The dependence of the NC cross section da/dQ"^, shown for the new e^p 
(solid points) and previously published 94 — 97 e~^p (open points) data. The error band and 
full curve represent the Standard Model expectations determined from the HI PDF 2000 fit at 
y/s = 319 GeV and v^s = 301 GeV, respectively, (b) The ratios of the 94 - 97 and 99 - 00 
data to their corresponding Standard Model expectations, where the normalisation shifts as 
determined from the fit are applied to the data (see table 3). The error band shows the Standard 
Model uncertainty for y/s = 319 GeV by adding in quadrature the experimental uncertainty as 
derived from the fit and the model uncertainty (see section 5.3). In (a) and (b), the inner and 
outer error bars represent respectively the statistical and total errors. The luminosity uncertainty 
is not included in the error bars. 
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Figure 6: (a) The dependence of the CC cross section da/dQ"^, shown for the new e^p 
(solid points) and previously published 94 — 97 e~^p (open points) data. The error band and 
full curve represent the Standard Model expectations determined from the HI PDF 2000 fit at 
y/s = 319 GeV and y/s = 301 GeV, respectively, (b) The ratios of the 94 - 97 and 99 - 00 
data to their corresponding Standard Model expectations, where the normalisation shifts as 
determined from the fit are applied to the data (see table 3). The error band shows the Standard 
Model uncertainty for y/s = 319 GeV. The error bars and band are defined as for fig. 5. 
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Figure 7: The dependences of the NC (circles) and CC (squares) cross sections da/dQ'^, 
shown for the combined 94 — 00 e^p (solid points) and 98 — 99 e^p (open points) data. The 
results are compared with the corresponding Standard Model expectations determined from the 
HI PDF 2000 fit. The error bars and bands are defined as for fig. 5. 



33 



Neutral Current 



a 3000 
2500 

■a 

2000 



1500 



1000 



500 











(a) 


Q^> 1000 GeV^ y<0.9 




* \ \ 1 j^^r» QQ nn a/c — 1 Q Ci(^\! 






o HI e'^p 94-97 Vs = 301 Ge\r 






+ -■ HlPDF2000Vs = 319GeV- 


- 




A — HI PDF 2000 Vs = 301 GeV- 




1 













10 



10 



J3 12 

Oh 



Neutral Current 



-a 



-a 



(b) 



10 - 



Q^> 10000 GeV^ y<0.9 
HI e'^p 99-00 Vs = 319 GeVj 
HI e*p 94-97 Vs = 301 GeV 
■ HI PDF 2000 Vs = 319 GeV - 
HI PDF 2000 Vs = 301 GeV^ 




Figure 8: The x dependence of the NC cross section da/dx for (a) Q'^ > 1 000 GeV^ and (b) 
> 10 000 GeV^, shown for the new e^p (solid points) and previously published 9A — 97 e+p 
(open points) data. The error bands and full curves represent the corresponding Standard Model 
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band and full curve represent the corresponding Standard Model expectations determined from 
the HI PDF 2000 m at ^/s = 319 GeV and ^ = 301 GeV, respectively The error bars and 
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Figure 10: The NC reduced cross section a^cix^ Q"^), shown for the new e^p (solid points) 
andhigh-y analysis of the 98 — 99 e~p (open triangles) data. The full (dashed) curves, labelled 
as gncqcd (o'NC,FL=o), show the Standard Model expectations detennined from the HI PDF 
2000 fit by including (excluding) the Fl contribution in the reduced cross section. The error 
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Figure 11: The NC reduced cross section aNci^-, Q"^)^ shown for the new e^p (solid points) 
and previously published 94 — 97 e+p (open points) data. The results are compared with the 
corresponding Standard Model expectations determined from the HI PDF 2000 fit at y/s = 
319 GeV (error bands) and ^/s = 301 GeV (full curves), respectively. Also shown are data 
from HI measured at lower Q"^ (open squares), as well as from the fixed-target experiment 
BCDMS (open triangles). The BCDMS data are not used in the fit. The error bars and bands 
are defined as for fig. 5. 
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Figure 12: The CC reduced cross section accix, Q"^), shown for the new e^p (solid points) 
and previously published 94 — 97 e~^p (open points) data. The results are compared with the 
corresponding Standard Model expectations determined from the HI PDF 2000 fit at a/s = 
319 GeV (fall curves) and y/s — 301 GeV (dashed curves), respectively. The error bars are 
defined as for fig. 5. 
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(solid points). The results are compared with the corresponding Standard Model expectation 
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Figure 14: Parton distributions (a) xU, (b) xU, (c) xD, (d) xD and (e) xg as determined from the 
HI PDF 2000 fit to HI data only. The distributions are shown at the initial scale Ql = 4: GeV^. 
The inner error band represents the experimental uncertainty as determined from the fit. The 
outer error band shows the total uncertainty by adding in quadrature the experimental and model 
uncertainties (see text). The valence quark distributions xuy (a) and xd^ (c) are also shown. For 
comparison, the parton distributions from the fit to HI and BCDMS data are shown as the full 
curves. 
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Figure 15: Parton distributions (a, b) xU and xU, (c, d) xD and xD, and (e, f) xg as determined 
from the HI PDF 2000 fit to HI data only. The distributions are shown at = 10 GeV^ (a, c, 
e) and at = i qoo GeV^ (b, d, f). The error bands are defined as for fig. 14. For comparison, 
recent results from the MRST and CTEQ groups are also shown. 
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Figure 16: The proton structure function F2 shown for the combined 94 — 00 e~^p (solid points) 
and previously published low (open circles) data. The results are compared with the corre- 
sponding Standard Model expectation determined from the HI PDF 2000 fit (error bands). The 
dashed curves show the backward extrapolation of the fit to < Also shown are the 

F2 data from BCDMS and NMC, which are not used in the fit. The error bars and bands are 
defined as for fig. 5. 
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Figure 17: Determination of Fl shown for the high y e~p and e^p data at fixed y = 0.75 as a 
function of (lower scale), or equivalently x (upper scale). The inner error bar represents the 
statistical error, the intermediate error bar shows the statistical and systematic errors added in 
quadrature and the outer error bar also includes the uncertainty arising from the extrapolation 
of F2. The error band, defined as for fig. 5, shows the expectation for Fl and its uncertainty, 
determined from the HI PDF 2000 fit. The full curve shows the expectation for Fl determined 
from the HI Low-y fit. The upper and lower dashed curves are the maximum and minimum 
allowed values for Fl. 



42 



u 



= 1500 GeV^ 



319 GeV 



0.75 - 



0.5 



0.25 - 




0.2 



0.1 



= 5000 GeV^ 



:: HI PDF 2000: 

+ 

e p 

e"p 



HI 

HI PDF 2000 





Q^= 12000 GeV^ 



• e""? HI 94-00 
o e"p HI 98-99 



(a) 




(b) 




10 



10 



10 




0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



Figure 18: The measured NC reduced cross sections d-^fj{x, Q^) (a), structure functions xF^ 
(b) and xF^^ (c), shown for three different values. The resuhs are compared with the cor- 
responding Standard Model expectations determined from the HI PDF 2000 fit. In (d), the 
averaged structure function xF]^ for a value of 1 500 GcV^ is compared with the expecta- 
tion determined from the same ht. The error bars are defined as for fig. 5. The normalisation 
uncertainties of the e~p and e^p data sets are included in the systematic errors. 
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Figure 19: The quark distributions xu and xd obtained from the local extraction method (solid 
points) in comparison with the expectations from the HI PDF 2000 fit (error bands), MRST 
(dashed curves) and CTEQ (dotted curves). 
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Table 6: The NC e^p cross section da^c/dQ^ for y < 0.9 after correction (kcor) according 
to the Standard Model expectation determined from the HI PDF 2000 fit for the kinematic cut 
E'^ > 6 GeV for < 890 GeV^. The statistical (Sgtat), uncorrelated systematic (Sunc), corre- 
lated systematic (5cor) and total (5tot) errors are provided. In addition the correlated systematic 
error contributions from a positive variation of one standard deviation of the electron energy 
error (S^^), of the polar electron angle error (6^^^), of the hadronic energy error (5^^^), of the 
error due to noise subtraction (S^Jl) and of the error due to background subtraction (5^^) are 
given. The normalisation uncertainty of 1.5% is not included in the errors. 
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5.1 


3.1 


9.4 


0.6 


2.7 


1.4 


0.0 


-12.0 


8000 


0.347 • 10-3 


1.04 


11.9 


9.4 


5.4 


16.1 


0.3 


5.1 


1.8 


0.0 


-11.1 


15000 


0.192 - 10-' 


1.06 


21.7 


16.9 


6.5 


28.2 


0.2 


6.1 


0.9 


0.0 


-15.7 



Table 7: The CC e+p cross section dacc/dQ^ fory < 0.9 after correction (kcor) according to 
the Standard Model expectation determined from the HI PDF 2000 fit for the kinematic cuts 
0.03 < y < 0.85 and Ptj,, > 12 GcV. The statistical (Sstat), uncorrelated systematic (5unc), 
correlated systematic (5cor) and total (5tot) errors are also given. In addition the correlated 
systematic error contributions from a positive variation of one standard deviation of the error 
due to the cuts against photoproduction (S^^l), of the hadronic energy error (5^^,.), of the error 
due to noise subtraction (S^J^) and of the error due to background subtraction (S^^) are given. 
The normalisation uncertainty of 1.5% is not included in the errors. The last column gives the 
correction for QED radiative effects (Aq^^). 



X 


daNc/dx (pb) 
g2 > 1000GeV2 


^stat 


^unc 


^cor 


5tot 


^cor 


^cor 


^cor 


^cor 


^cor 




(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 




y < 0.9 




















0.013 


0.138 • lO'' 


5.9 


4.2 


3.7 


8.1 


-0.8 


-1.0 


-1.5 


-0.8 


-3.1 


0.020 


0.249 - 10^ 


3.0 


2.3 


1.4 


4.1 


0.3 


-1.2 


0.4 


0.3 


-0.6 


0.032 


0.215 • 10"^ 


2.6 


2.3 


1.2 


3.6 


0.6 


-0.7 


0.7 


0.4 


-0.1 


0.050 


0.147- 10"^ 


2.5 


2.0 


1.0 


3.3 


0.3 


-0.7 


0.5 


0.4 


-0.1 


0.080 


0.951 - 10^ 


2.4 


2.1 


1.1 


3.4 


0.5 


-0.6 


0.5 


0.6 


0.0 


0.130 


0.566 - 10^ 


2.8 


2.2 


1.1 


3.7 


0.9 


-0.5 


0.3 


0.4 


0.0 


0.180 


0.372 • 10^ 


3.0 


2.6 


1.6 


4.3 


1.3 


-0.8 


-0.3 


0.5 


0.0 


0.250 


0.212 - 10'^ 


3.3 


3.4 


2.1 


5.2 


2.0 


-0.2 


-0.6 


-0.3 


0.0 


0.400 


0.646 - 102 


4.7 


5.7 


5.3 


9.1 


3.7 


-0.4 


-1.9 


-3.5 


0.0 


0.650 


0.650 - lOi 


9.8 


13.3 


11.6 


20.2 


7.4 


0.8 


-4.0 


-8.1 


0.0 



Table 8: The NC e^p cross section daNc/dx measured for y < 0.9 and > 1 000 GeV^. 
The statistical (5stat), uncorrelated systematic (Sunc), correlated systematic (Scar) and total (5tot) 
errors are provided. In addition the correlated systematic error contributions from a positive 
variation of one standard deviation of the electron energy error (S^^), of the polar electron 
angle error (5^^^), of the hadronic energy error (6'^^^.), of the error due to noise subtraction (5^) 
and of the error due to background subtraction (5^^r) given. The normalisation uncertainty 
of 1.5% is not included in the errors. 
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X 


daNc/dx (pb) 
Q2 > lOOOOGeV^ 
y < 0.9 


^stat 
(%) 


^unc 
(%) 


(%) 


^tot 

(%) 


"cor 
(%) 


"cor 
(%) 


6^^ 

^cor 
(%) 


"cor 
(%) 


6^^ 

^cor 
(%) 


0.130 


0.544 • 10^ 


29.7 


8.9 


1.9 


31.0 


0.1 


-1.2 


-1.1 


-0.6 


-0.8 


0.180 


0.716 • 10^ 


18.9 


4.2 


2.2 


19.5 


-1.1 


-1.9 


0.4 


0.1 


-0.5 


0.250 


0.411 • 10^ 


19.7 


6.3 


1.4 


20.7 


1.1 


-0.5 


0.9 


0.3 


-0.1 


0.400 


0.162 • 10^ 


24.5 


11.5 


2.0 


27.1 


1.9 


0.6 


0.5 


0.2 


-0.3 


0.650 


0.038 • 10^ 


31.6 


30.5 


8.4 


44.8 


6.3 


5.1 


-2.2 


-0.8 


0.0 



Table 9: The NC e^p cross section daNc/(ix measured for y < 0.9 and Q"^ > 10 000 GeV^. 
The statistical (dgtat), uncorrelated systematic (5unc), correlated systematic (Scor) and total 0tot) 
errors are provided. In addition the correlated systematic error contributions from a positive 
variation of one standard deviation of the electron energy error of the polar electron 

angle error (^^or)' hadronic energy error (5^o^.), of the error due to noise subtraction (5^J,.) 
and of the error due to background subtraction (S^J^) are given. The normalisation uncertainty 
of 1.5% is not included in the errors. 



X 


d(Tcc/dx (pb) 
g2 > 1 000 GeV^ 
y < 0.9 


k 

"'Cor 


Sstat 
(%) 


^unc 
(%) 


Scor 
(%) 


Stot 

(%) 


"cor 
(%) 


"cor 
(%) 


^7V+ 
"cor 

(%) 


"cor 
(%) 


0.032 


0.159 • 10^ 


1.05 


6.9 


4.0 


2.3 


8.3 


2.1 


-0.5 


0.6 


-0.4 


0.080 


0.946 • 10^ 


1.02 


4.6 


3.4 


1.2 


5.9 


0.9 


0.5 


0.7 


0.0 


0.130 


0.623 ■ 10^ 


1.01 


5.2 


3.7 


2.0 


6.7 


0.3 


1.2 


1.5 


0.0 


0.250 


0.194-102 


1.00 


7.2 


5.8 


2.2 


9.5 


0.1 


2.2 


0.1 


0.0 


0.400 


0.451 • 10^ 


1.06 


16.6 


11.3 


6.0 


21.0 


0.0 


4.9 


-3.6 


0.0 


0.650 


0.469 • 10° 


1.26 


70.6 


24.5 


21.3 


77.7 


0.0 


9.2 


-19.2 


0.0 



Table 10: The CCe^p cross section dacc/dx fory < 0.9 andQ^ > 1 000 GeV^ after correction 
(kcor) according to the Standard Model expectation determined from the HI PDF 2000 fit for 
the kinematic cuts 0.03 < y < 0.85 and Pt^h > 12 GeV. The statistical (Sgtat), uncorrelated 
systematic (5unc), correlated systematic (Scor) and total (Stot) errors are also given. In addition 
the correlated systematic error contributions from a positive variation of one standard deviation 
of the error due to the cuts against photoproduction (S^orX of the hadronic energy error (6^^^.), 
of the error due to noise subtraction (5^^) and of the error due to background subtraction (S^ot) 
are given. The normalisation uncertainty of 1.5% is not included in the errors. 
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Table 11: The NC e+p reduced cross section aNc{x, Q'^), shown with statistical (6stat), systematic (dgys) and total (Stot) errors. 
Also shown are the total uncorrelated systematic (5unc) errors and two of its contributions: the electron energy error (S^^c) 
hadronic energy error (S^nJ- effect of the other uncorrelated systematic errors is included in (Sunc)- The table also provides 
the correlated systematic error (5f.(„.) and its contributions from a positive variation of one standard deviation of the error on the 
electron energy (6fJ^.) and polar angle (5cor)' hadronic energy error (^cor)' error due to noise subtraction (5^^) and 

background subtraction (5^^^) and of the error due to charge symmetry background subtraction in the high-y analysis (S^g^)- The 
normalisation uncertainty of 1.5% is not included in the errors. The NC structure function term scaled by the helicity factor 
4'Nc/Y+ is given as well as the electromagnetic structure function F2 with the corrections Ap.^, Ap.^ and Ap^ as defined in eg. 18. 
For Q'^ < 2 000 GeV^, the extraction of F2 is restricted to the region ofy< 0.6. The table continues on the next 2 pages. 
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0.65000 


0.023 


0.0166 


30.2 


19.4 


35.9 


14.4 


8.4 


5.6 


13.0 


8.5 


-2.7 


-3.9 


-8.7 


0.0 




0.0166 


0.0160 


4.4 


-0.4 


0.0 


2000 


0.03200 


0.615 


0.642 


6.4 


4.2 


7.7 


3.9 


0.5 


1.5 


1.6 


0.6 


— 0.5 


1.3 


0.5 


0.0 




0.644 


0.671 


3.8 


— 5 . 3 


— 2.5 


2000 


0.05000 


0.394 


0.472 


7.4 


4.0 


8.4 


3.7 


0.9 


0.9 


1.7 


1.2 


-0.6 


0.8 


0.6 


0.0 




0.473 


0.481 


3.8 


-4.7 


-0.7 


2000 


0.08000 


0.246 


0.405 


6.7 


3.5 


7.6 


3.2 


0.4 


0.9 


1.5 


-0.4 


-1.0 


0.8 


0.7 


0.0 




0.406 


0.408 


3.5 


-3.8 


-0.2 


2000 


0.13000 


0.151 


0.368 


7.4 


4.0 


8.4 


3.6 


1.5 


0.4 


1.7 


1.6 


-0.8 


0.1 


0.0 


0.0 




0.369 


0.367 


3.3 


-2.8 


0.0 


2000 


0.18000 


0.109 


0.342 


8.3 


3.9 


9.2 


3.7 


0.3 


0.6 


1.0 


0.3 


0.2 


-0.8 


0.5 


0.0 




0.343 


0.340 


3.1 


-2.2 


0.0 


2000 


0.25000 


0.079 


0.233 


9.0 


4.9 


10.2 


4.4 


2.1 


1.0 


2.3 


2.2 


— 0.3 


— 0.6 


-0.3 


0.0 




0.233 


0.230 


3.0 


-1.7 


0.0 


2000 


0.40000 


0.049 


0.108 


12.4 


10.1 


16.0 


7.5 


5.3 


2.6 


6.8 


5.3 


-0.9 


-2.4 


-3.3 


0.0 




0.108 


0.106 


2.8 


-1.1 


0.0 


2000 


0.65000 


0.030 


0.0294 


23.6 


20.5 


31.3 


14.7 


7.4 


6.9 


14.2 


7.4 


-1.4 


-3.9 


-11.4 


0.0 




0.0294 


0.0287 


3.2 


-0.7 


0.0 


3000 


0.05000 


0.591 


0.551 


5.6 


3.9 


6.8 


3.6 


0.4 


1.4 


1.5 


0.0 


-0.7 


1.2 


0.6 


-0.1 




0.553 


0.585 


6.3 


-9.9 


-1.7 


3000 


0.08000 


0.369 


0.439 


6.2 


4.0 


7.4 


3.6 


0.5 


1.1 


1.8 


1.2 


—0.5 


0.9 


0.8 


0.0 




0.441 


0.452 


6.0 


-8.2 


-0.4 


3000 


0.13000 


0.227 


0.340 


7.3 


3.8 


8.2 


3.7 


1.0 


0.5 


0.9 


0.8 


-0.4 


-0.4 


0.1 


0.0 




0.341 


0.343 


5.7 


-6.1 


-0.1 


3000 


0.18000 


0.164 


0.313 


7.8 


5.9 


9.8 


5.0 


3.3 


1.0 


3.2 


3.1 


-0.3 


-0.5 


0.1 


0.0 




0.314 


0.312 


5.3 


-4.7 


0.0 


3000 


0.25000 


0.118 


0.195 


9.7 


4.7 


10.8 


4.0 


1.1 


0.7 


2.4 


2.3 


0.4 


-0.7 


0.0 


0.0 


_ 


0.195 


0.192 


5.1 


-3.6 


0.0 


3000 


0.40000 


0.074 


0.108 


11.5 


6.4 


13.1 


5.5 


2.8 


1.4 


3.3 


2.9 


0.3 


-1.2 


-1.1 


0.0 


- 


0.108 


0.106 


4.6 


-2.3 


0.0 


3000 


0.65000 


0.045 


0.0118 


28.9 


16.8 


33.4 


12.4 


6.7 


5.4 


11.4 


6.6 


1.2 


-3.5 


-8.6 


0.0 


- 


0.0118 


0.0113 


5.7 


-1.4 


0.0 


5000 


0.08000 


0.615 


0.422 


7.0 


4.3 


8.2 


4.1 


0.9 


1.6 


1.5 


0.1 


-0.6 


1.3 


0.5 


-0.1 




0.425 


0.475 


11.0 


-20.1 


-1.4 


5000 


0.13000 


0.379 


0.339 


8.3 


4.3 


9.3 


4.3 


1.2 


1.0 


0.9 


0.6 


-0.2 


0.6 


0.2 


0.0 




0.341 


0.359 


10.4 


-15.2 


-0.3 


5000 


0.18000 


0.273 


0.263 


9.4 


4.7 


10.5 


4.7 


2.2 


0.3 


0.5 


-0.5 


-0.1 


-0.1 


0.0 


0.0 




0.263 


0.269 


9.9 


-11.8 


-0.1 


5000 


0.25000 


0.197 


0.301 


9.6 


9.4 


13.4 


9.0 


7.5 


1.3 


2.9 


2.8 


-0.5 


-0.6 


0.3 


0.0 




0.302 


0.300 


9.4 


-8.8 


0.0 


5000 


0.40000 


0.123 


0.130 


12.3 


6.1 


13.7 


5.9 


2.5 


0.3 


1.2 


1.1 


0.2 


0.3 


0.5 


0.0 




0.130 


0.127 


8.3 


-5.5 


0.0 


5000 


0.65000 


0.076 


0.00760 


40.9 


25.0 


47.9 


21.3 


17.1 


7.9 


13.2 


8.9 


2.4 


-6.3 


-7.1 


0.0 




0.00761 


0.00721 


8.9 


-3.4 


0.0 


8000 


0.13000 


0.606 


0.344 


11.0 


5.5 


12.3 


5.4 


1.8 


1.2 


1.0 


0.6 


0.3 


0.7 


0.1 


-0.1 




0.347 


0.419 


17.0 


-33.1 


-0.9 


8000 


0.18000 


0.438 


0.361 


10.7 


6.5 


12.5 


6.3 


3.1 


1.7 


1.6 


0.9 


-1.0 


0.9 


0.1 


0.0 




0.364 


0.406 


15.8 


-26.1 


-0.3 


8000 


0.25000 


0.315 


0.224 


13.1 


7.1 


14.9 


6.7 


3.9 


0.0 


2.5 


2.4 


0.4 


0.2 


0.7 


0.0 




0.224 


0.234 


15.2 


-19.4 


-0.1 


8000 


0.40000 


0.197 


0.0982 


18.0 


10.3 


20.7 


9.8 


7.2 


1.2 


3.2 


3.2 


-0.3 


-0.3 


0.0 


0.0 




0.0984 


0.0971 


13.4 


-12.0 


0.0 


8000 


0.65000 


0.121 


0.0228 


28.9 


28.7 


40.7 


26.6 


22.7 


6.3 


11.0 


8.2 


4.1 


-4.7 


-4.0 


0.0 




0.0228 


0.0213 


14.5 


-7.3 


0.0 


12000 


0.18000 


0.656 


0.233 


19.5 


5.0 


20.1 


4.6 


2.8 


0.4 


2.0 


-0.9 


-1.8 


0.4 


0.2 


-0.3 




0.236 


0.320 


22.7 


-48.1 


-0.9 


12000 


0.25000 


0.473 


0.117 


23.6 


7.3 


24.7 


7.1 


6.1 


1.5 


1.5 


1.0 


0.6 


0.9 


0.3 


0.0 




0.117 


0.139 


21.1 


-36.7 


-0.3 


12000 


0.40000 


0.295 


0.0424 


35.4 


11.8 


37.3 


11.6 


11.1 


0.2 


2.0 


1.8 


0.7 


-0.3 


0.1 


0.0 




0.0425 


0.0441 


19.1 


-22.6 


-0.1 


12000 


0.65000 


0.182 


0.0179 


40.8 


28.9 


50.0 


27.4 


26.2 


4.8 


9.0 


6.4 


5.1 


-3.6 


-1.1 


0.0 




0.0180 


0.0170 


19.1 


-13.5 


0.0 


20000 


0.25000 


0.788 


0.106 


38.3 


5.2 


38.7 


4.8 


2.9 


0.6 


2.0 


1.1 


-1.5 


0.9 


0.1 


-0.4 




0.109 


0.192 


29.5 


-71.8 


-1.0 


20000 


0.40000 


0.492 


0.0381 


50.0 


11.7 


51.4 


11.5 


11.0 


1.6 


2.4 


2.1 


0.7 


0.9 


0.3 


0.0 




0.0385 


0.0490 


26.7 


-48.0 


-0.2 


20000 


0.65000 


0.303 


0.0110 


70.7 


38.7 


80.6 


38.0 


37.5 


1.0 


7.7 


5.9 


4.9 


-0.7 


-0.3 


0.0 




0.0110 


0.0113 


26.2 


-28.6 


0.0 


30000 


0.40000 


0.738 


0.164 


46.1 


13.7 


48.1 


13.4 


12.6 


0.9 


2.4 


-0.6 


-2.2 


-0.4 


0.2 


-0.9 




0.168 


0.319 


32.4 


-79.1 


-0.5 



(GeV2) 


X 


y 


(pb/GeV2) 


4>cc 


^stat 
(%) 


^sys 

(%) 


Stot 

(%) 


^unc 

(%) 


(5'^ 

unc 

(%) 


(%) 


^cor 

(%) 


cor 

(%) 


^cor 

(%) 


^cor 

(%) 


^cc 

(%) 


300 

300 
300 


0.0130 

0.0320 
0.0800 


0.227 

0.092 
0.037 


0.703 • lO'J 
0.283 • 10° 
0.585-10-1 


1.184 
1.171 
0.606 


20.3 
13.7 
19.0 


11.8 
5.9 
7.4 


23.5 
14.9 
20.4 


8.2 
4.6 
6.0 


1.5 
1.3 
2.9 


8.6 
3.8 
4.5 


6.7 
2.6 
1.0 


-2.2 
-1.4 
-2.4 


-0.3 
1.0 
1.6 


-4.9 
-2.1 
-3.2 


0.3 
0.4 
5.2 


500 
500 
500 
500 


0.0130 
0.0320 
0.0800 
0.1300 


0.379 
0.154 
0.062 
0.038 


0.570 • 10° 
0.189-10° 
0.465 - 10-1 
0.194-10-1 


1.018 
0.829 
0.511 
0.346 


14.6 
12.1 
13.4 
25.1 


8.5 

4.8 
4.4 
7.2 


16.9 
13.0 
14.1 
26.1 


6.5 
3.9 
4.0 
6.4 


2.3 
1.2 
0.4 
2.2 


5.4 
2.8 
2.0 
2.9 


4.9 
1.8 
0.4 
0.2 


-1.4 
-1.1 
-0.7 
-1.6 


0.4 
0.8 
1.9 

-2.2 


-1.8 
-1.6 
-0.1 
-0.9 


-4.4 
-0.7 
-0.7 
-3.5 


1000 
1000 
1000 
1000 


0.0320 
0.0800 
0.1300 
0.2500 


0.308 
0.123 
0.076 
0.039 


0.121 - 10° 
0.406 - 10-1 
0.162-10-1 
0.794-10-2 


0.609 
0.512 
0.332 
0.313 


10.5 
10.2 
16.5 
23.5 


4.7 
3.5 
6.1 
16.8 


11.5 
10.8 
17.6 

28.9 


3.8 
3.3 
5.9 
16.1 


1.3 
0.6 
0.5 
2.0 


2.5 
1.2 
0.9 

4.9 


1.9 
0.6 
0.3 
0.0 


-1.3 
-0.7 
0.5 
2.4 


0.4 
0.8 
0.7 
-4.2 


-0.8 
0.0 
0.0 
0.0 


-3.1 
-0.4 
-2.5 
-4.9 


2000 
2000 
2000 
2000 


0.0320 
0.0800 
0.1300 
0.2500 


0.615 
0.246 
0.152 
0.079 


0.762 - 10-1 
0.228- 10-1 
0.168-10-1 
0.337-10-2 


0.495 
0.370 
0.442 
0.171 


9.8 
9.9 
11.1 
18.4 


4.5 
3.8 
5.7 
7.3 


10.8 
10.6 
12.5 
19.8 


4.0 
3.6 
5.2 
6.9 


0.2 
0.6 
2.0 
0.3 


2.3 
1.2 
2.6 
2.0 


2.1 
0.8 
0.1 
0.0 


0.1 
-0.8 
-1.3 
-1.4 


0.8 
0.4 
2.2 
-1.5 


-0.2 
0.0 
0.0 
0.0 


-5.6 
-3.1 
-5.0 
-10.0 


3000 
3000 
3000 
3000 


0.0800 
0.1300 
0.2500 
0.4000 


0.369 
0.227 
0.118 
0.074 


0.201 - 10-1 
0.107-10-1 
0.251 - 10-2 
0.531 - 10-3 


0.407 
0.354 
0.159 
0.054 


8.7 
10.7 
16.3 
37.8 


4.3 
4.5 
6.1 
17.1 


9.7 
11.6 
17.4 
41.5 


4.2 
3.8 
6.0 
14.2 


1.2 
1.1 
0.4 
3.9 


1.4 
2.1 
1.1 
9.6 


0.9 
0.2 
0.1 
0.0 


0.9 
1.2 
1.0 
1.8 


0.2 
1.7 
-0.4 
-9.4 


0.0 
0.0 
0.0 
0.0 


-7.5 
-3.9 
-5.7 
-12.6 


5000 
5000 
5000 
5000 


0.0800 
0.1300 
0.2500 
0.4000 


0.615 
0.379 
0.197 
0.123 


0.842 - 10-2 
0.530-10-2 
0.192- 10-2 
0.261 • 10-3 


0.250 
0.256 
0.179 
0.039 


13.8 
12.3 
14.2 
33.3 


8.1 
5.6 
5.7 
10.4 


16.0 
13.5 
15.3 
34.9 


5.9 

5.1 
5.3 
9.5 


1.8 

2.5 
3.0 
4.0 


5.6 
2.3 
1.7 

4.6 


1.4 
0.4 
0.1 
0.0 


3.8 
2.2 
1.4 
4.5 


3.8 
0.5 
0.9 

-0.5 


0.0 
0.0 
0.0 
0.0 


-13.0 

-13.8 
-9.8 
-4.6 


8000 
8000 
8000 


0.1300 
0.2500 
0.4000 


0.606 
0.315 
0.197 


0.178-10-^ 
0.903-10-3 
0.152-10-3 


0.137 
0.134 
0.036 


20.3 
17.3 
40.8 


12.6 
11.0 
26.0 


23.9 
20.5 
48.4 


11.4 
9.5 
25.2 


5.1 
8.1 
14.5 


5.4 
5.8 
6.6 


0.6 
0.0 
0.0 


4.7 
5.3 
6.6 


2.7 
2.2 
0.2 


0.0 
0.0 
0.0 


-13.0 
-8.6 

-14.8 


15000 
15000 


0.2500 
0.4000 


0.591 
0.369 


0.126-10-3 
0.117- 10-3 


0.041 
0.061 


37.8 
28.8 


16.4 
23.2 


41.2 
37.0 


15.0 
22.0 


13.6 
13.4 


6.8 
7.6 


0.3 
0.0 


6.5 
7.4 


1.7 
1.3 


0.0 
0.0 


-16.0 
-14.7 



Table 12: The CC e^p double differential cross section d^acc/dxdQ^ and the structure function term (f)cc, shown 
with statistical (SgtatX systematic (Sgys) and total (5tot) errors. Also shown are the total uncorrelated systematic error 
(dune) its contribution from the hadronic energy error (Sll^^J. The effect of the other uncorrelated systematic errors 
is included in (5y,nc)- In addition the correlated systematic error (Scor) and its contributions from a positive variation 
of one standard deviation of the error due to the cuts against photoproduction (5^^), of the hadronic energy error 
(5^^,.), of the error due to noise subtraction (5^^.) and of the error due to background subtraction (S^J,.) are given. The 
normalisation uncertainty of 1.5% is not included in the errors. The last column gives the correction for QED radiative 
effects (Agg^). 



r)2 
Q 


X 




Otot 


Ostat 






gh 


dear 


SE + 




"cor 


gN + 


gB + 


"cor 


(GeV^) 






{%) 


{%) 


{%) 


{%) 


(%) 


(%) 


{%) 


{%) 


(%) 


(%) 


{%) 


(%) 


100 


0.00130 


1.368 


7.0 


4.0 


5.5 


1.5 


0.2 


1.2 


-0.7 


-1.0 


0.5 


0.4 


- 


-0.9 


100 


0.00200 


1.293 


5.3 


3.4 


3.3 


0.5 


0.3 


2.4 


0.5 


-1.8 


0.9 


0.6 


-1.1 




120 


0.00160 


1.342 


6.6 


4.2 


4.8 


0.1 


0.2 


1.7 


-0.3 


-1.5 


0.6 


0.3 




-0.7 


120 


0.00200 


1.325 


5.0 


3.3 


3.2 


0.5 


0.2 


2.1 


-0.4 


-1.4 


-0.3 


0.3 


-1.4 




120 


0.00320 


1.198 


4.8 


3.1 


3.3 


1.1 


0.4 


1.9 


0.8 


-1.7 


0.5 


0.5 


-0.3 




150 


0.00200 


1.339 


6.7 


4.4 


4.8 


0.9 


0.2 


1.2 


0.8 


-0.8 


0.5 


0.4 




-0.6 


200 


0.00260 


1.188 


7.1 


4.9 


4.8 


0.1 


0.4 


1.9 


0.4 


-1.7 


0.6 


0.4 




-0.6 


250 


0.00330 


1.126 


7.9 


5.7 


4.9 


0.6 


0.2 


2.1 


1.0 


-2.0 


0.5 


0.3 




-0.7 


:«)() 


11. ()():!<)() 


l.()(i,^ 


S.ll 


(i.l 


4.!) 


0.2 


11.2 


1.4 


0.4 


i . :i 


11.4 


11.3 




-0.7 


400 


u.()()5:i() 


l.iOl 




(>.4 


5.1 


0.7 


u.:i 


1.5 


-0.6 


-1.4 


o.:i 


0.3 




-0.6 


500 


0.00660 


1.099 


8.7 


6.9 


5.2 


0.2 


0.3 


1.0 


-0.3 


-0.8 


0.5 


0.4 




-0.2 


650 


0.00850 


1.056 


9.9 


8.2 


5.6 


1.0 


0.0 


0.5 


-1.0 


-0.2 


-0.2 


0.2 




-0.4 


800 


0.01050 


0.938 


10.9 


9.2 


5.9 


0.4 


0.3 


1.1 


-0.6 


-1.1 


0.3 


0.0 




-0.1 



Table 13: The JVC e~p reduced cross section a^cix^Q'^) from the high-y analysis, shown with 
statistical (5stat) and total (5tot) errors. Also shown are the total uncorrelated systematic (5unc) 
error and two of its contributions: the electron energy error (S^^c) hadronic energy 

error (S'^^J. The effect of the other uncorrelated systematic errors is included in (Sunc)- 
addition the correlated systematic error (5a„.) and its contributions from a positive variation of 
one standard deviation of the electron energy error (5^^), of the polar electron angle error (5^^^), 
of the hadronic energy error (6^^^), of the error due to noise subtraction (5^^), of the error due to 
background subtraction (5^^) and of the error due to charge symmetry background subtraction 
(S^oj.) are given. The normalisation uncertainty of 1.8% is not included in the errors. All e~p 
data not previously reported in [3] are given, including the new high y data and three data points 
at = 100 and 120 GeV^ from the nominal analysis phase space. 



g2 
(GeV^) 


X 


y 


(pNC 


Fl 




^sys 


5tot 


e p data 


110 


0.00144 


0.75 


1.440 


0.298 


0.074 


0.133 


0.154 


175 


0.00230 


0.75 


1.346 


0.298 


0.077 


0.113 


0.139 


280 


0.00368 


0.75 


1.162 


0.390 


0.085 


0.099 


0.132 


450 


0.00591 


0.75 


1.164 


0.117 


0.097 


0.101 


0.140 


700 


0.00919 


0.75 


1.072 


0.042 


0.117 


0.098 


0.153 


e^p data 


110 


0.00144 


0.75 


1.518 


0.198 


0.032 


0.083 


0.092 


175 


0.00230 


0.75 


1.426 


0.164 


0.038 


0.064 


0.076 


280 


0.00368 


0.75 


1.292 


0.171 


0.041 


0.057 


0.072 


450 


0.00591 


0.75 


1.163 


0.133 


0.045 


0.052 


0.070 


700 


0.00919 


0.75 


1.037 


0.096 


0.053 


0.062 


0.082 



Table 14: The NC structure function term (f)]sfc{x, Q"^) and the structure function F^, shown with 
its statistical (SgtatX systematic (days) and total (5tot) absolute error. The total error includes 
a contribution arising from the model uncertainties in the calculated F2. These are obtained 
by varying the assumptions of the HI Low y QCD fit as listed in table 5. The luminosity 
uncertainties of the e^p and e~p data sets are included in the systematic error. 
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(GeV2) 


X 




^stat 


^sys 


^tot 


1500 


0.020 


0.052 


0.036 


0.025 


0.044 


1500 


0.032 


0.074 


0.032 


0.026 


0.042 


1500 


0.050 


0.076 


0.039 


0.028 


0.048 


1500 


0.080 


0.067 


0.050 


0.035 


0.061 


5000 


0.050 


0.088 


0.037 


0.024 


0.044 


5000 


0.080 


0.150 


0.031 


0.020 


0.037 


5000 


0.130 


0.160 


0.036 


0.023 


0.043 


5000 


0.180 


0.099 


0.041 


0.025 


0.048 


5000 


0.250 


0.089 


0.049 


0.039 


0.062 


5000 


0.400 


0.027 


0.045 


0.034 


0.057 


5000 


0.650 


-0.008 


0.015 


0.011 


0.019 


12000 


0.180 


0.149 


0.077 


0.021 


0.080 


12000 


0.250 


0.113 


0.053 


0.017 


0.056 


12000 


0.400 


0.035 


0.038 


0.019 


0.043 


12000 


0.650 


-0.007 


0.015 


0.011 


0.018 


(GeV2) 


X 


xF2' 


^stat 


^sys 


Stot 


1500 


0.026 


0.59 


0.22 


0.17 


0.28 


1500 


0.050 


0.38 


0.13 


0.09 


0.16 


1500 


0.080 


0.57 


0.12 


0.08 


0.14 


1500 


0.130 


0.61 


0.14 


0.09 


0.16 


1500 


0.180 


0.37 


0.12 


0.06 


0.13 


1500 


0.250 


0.29 


0.11 


0.05 


0.12 


1500 


0.400 


0.09 


0.08 


0.05 


0.09 


1500 


0.650 


-0.02 


0.03 


0.02 


0.04 



Table 15: The upper part of the table shows the generalised structure function xF^ with statisti- 
cal (5 Stat), systematic (6 sys) and total (Stot) absolute errors. The luminosity uncertainties of the 
e^p and e~p data are included in the systematic error. The lower part of the table shows the 
structure function xF^^ obtained by averaging over different values and transforming to a 
g2 value at 1500 GeVl 
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